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ABSTRACT 

Context. Mass loss from red giants in old globular clusters affects the horizontal branch (HB) morphology and post-HB stellar 
evolution including the production of ultraviolet-bright stars, dredge up of nucleosynthesis products and replenishment of the intra- 
cluster medium. Studies of mass loss in globular clusters also allows one to investigate the metallicity dependence of the mass loss 
from cool, low-mass stars down to very low metallicities. 

Aims. We present an analysis of new VLT/UVES spectra of 47 red giants in the Galactic globular clusters 47 Tuc (NGC 104), NGC 
362, CO Cen (NGC 5139), NGC 6388, M54 (NGC 6715) and M15 (NGC 7078). The spectra cover the wavelength region 6100-9900 A 
at a resolving power ofR = 110, 000. Some of these stars are known to exhibit mid-infrared excess emission indicative of circumstellar 
dust. Our aim is to detect signatures of mass loss, identify the mechanism(s) responsible for such outflows, and measure the mass-loss 
rates. 

Methods. We determine for each star its effective temperature, luminosity, radius and escape velocity. We analyse the Ha and near- 
infrared calcium triplet lines for evidence of outflows, pulsation and chromospheric activity, and present a simple model for estimating 
mass-loss rates from the Ho- line profile. We compare our results with a variety of other, independent methods. 
Results. We argue that a chromosphere persists in Galactic globular cluster giants and controls the mass-loss rate to late-K/early-M 
spectral types, where pulsation becomes strong enough to drive shock waves at luminosities above the RGB tip. This transition may 
be metallicity-dependent. We find mass-loss rates of ~10"^ to 10"^ Mq yr"', largely independent of metallicity. 

Key words. Stars: AGB and post-AGB ~ Stars: circumstellar matter - Stars: mass-loss - Stars: chromospheres - Infrared: stars - 
Galaxy: globular clusters: individual: NGC 104, NGC 362, NGC 5139, NGC 6388, NGC 6715, NGC 7078 



1. Introduction 

. Understanding the evolution and mass loss in red giant branch 
(RGB) and asymptotic giant branch (AGB) stars is of substan- 
tial importance in uncovering the history of metal enrichment of 
the inter-stellar medium (ISM) and subsequent stellar and plan- 
etary formation; as well as understanding the evolution of the 
stars themselves. By observing mass-losing giant branch stars 
in globular clusters (GCs), we can simultaneously observe a co- 
eval set of objects within a cluster with similar progenitor mass 
and identical metallicity; while we can use different clusters to 
observe the effects of differing ages and masses, but more partic- 
ularly metallicities, which can range from near-solar metallicity 
to as little as ~Q.5% of solar. 

Simultaneously, we also glean information about the glob- 
ular cluster environment itself - the mass lost from stars forms 
the intra-cluster medium (ICM). Studies of the ICM show it is 
removed from the cluster on a timescale much shorter than the 
time between Galactic Plane crossings (e.g. Boyer et al. 2006; 
van Loon et al. 2006a). This mass loss exacerbates cluster evapo- 
ration, possibly leading to their dispersal on a timescale of 10^ 
years (van Loon & McDonald 2007). 

Both the RGB and AGB are associated with episodes of 
strong mass loss, eventually leading to stellar death. Globular 
cluster stars are thought to lose around 0.2 M© over their time on 
the RGB (Rood 1973) with up to 0. 1 Mq being lost in the helium 
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flash at the RGB tip (Fusi-Pecci & Renzini 1975, and references 
therein). This is required to explain the horizontal branch (HB) 
morphology, which is dependent on the remaining mantle mass. 
Dupree et al. (2007) have now observationally confirmed mass 
loss three magnitudes below the RGB tip, thought to be driven by 
hydrodynamically- or acoustically-dominated chromospheres. 

Stars reaching the upper slopes of the AGB are even more 
extreme, losing mass at rates reaching well over 10"* M© yr"' 
(Wood et al. 1983, 1992; van Loon et al. 1999). This is facilitated 
by thermal pulses - periodically occurring temporary ignition 
of a helium shell-burning source inside the star - and shorter- 
timescale (~10^"^ days) radial pulsations at the surface which 
can lead to shock-wave emission in the extended stellar atmo- 
sphere (Fox & Wood 1985). In metal-rich AGB stars, dust form- 
ing in the outmost parts of these atmospheres is subject to radi- 
ation pressure, which forces it and the gas it collides with to be 
radially accelerated away from the star It is not yet clear if dust- 
driving is the dominating effect in RGB and metal-poor AGB 
winds, and it is possible that Alfven waves or pressure (P-mode) 
waves dominate here (Hartmann & MacGregor 1980; Pijpers & 
Habing 1989). Indeed, Judge & Stencel (1991, hereafter JS91) 
argue that the mass-loss rate is largely invariant with the driv- 
ing process and that most cool giants (RGB or AGB) do not 
have dust as the primary driver of their winds. More recently, 
Schroder & Cuntz (2005, hereafter SC05) have suggested an 
improved Reimers relation (Reimers 1975), by assuming this in- 
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variance and using simple chromospheric energy considerations, 
which is calibrated using the total RGB mass-loss of GC stars. 

Previous works have attempted to provide evidence for mass 
loss and calculate mass-loss rates by analysing optical spectral 
line profiles, either through asymmetries in the line cores (first 
performed by Deutsch (1956), for the field giant a Her), line 
emission (Cohen 1976) or modelling. Line emission has also 
been used to model a chromosphere with a large bulk motion 
(Dupree et al. 1984; Mauas et al. 2006). Semi-empirical methods 
of finding mass-loss rates have also been used to define relation- 
ships based on stellar parameters. Mass-loss rates can also be de- 
rived from dust masses calculated from IR excess emission (e.g. 
Origlia et al. 2002 - hereafter OFFR02; van Loon 2007), though 
this must assume a value for the ratio of the masses of dust and 
gas lost. Mira-type long-period variability is only seen in stars 
with [Fe/H] > -1 (Frogel & Whitelock 1998), so one would 
expect that, assuming pulsation assists in driving mass-loss, the 
mass-loss rate and/or mechanism would vary with metallicity. 

In this work, we use all of the above techniques to try to iden- 
tify mass loss and, where possible, estimate mass-loss rates from 
a sample of 47 such objects from a range of clusters. Ultimately, 
we aim to assess the role of dust, pulsation and chromospheric 
activity in driving mass loss, and its dependence on metallicity. 
Our observations are based on new VLT/UVES echelle spectra 
at very high spectral resolution. We present the spectra and liter- 
ature photometry in Section 2 of the paper. In Section 3, we use 
the Kurucz atlas9 models (Kurucz 1993) to estimate the effec- 
tive temperature, metallicity and surface gravity of the stars, then 
use the near-IR photometry to calculate other physical parame- 
ters (luminosity, radius and escape velocity). A full abundance 
analysis has not been performed, but will become the subject of 
a subsequent paper In Section 4, we analyse the Ha and near- 
IR Ca II triplet lines and from them provide evidence of outflow 
from the stellar surface. We present a simple model for estimat- 
ing mass loss in the stellar wind from the Ha profile in Section 
5, and use it to calculate mass-loss rates and wind velocities for 
our sample of stars. This is followed by a discussion in Section 6, 
in which we also estimate mass-loss rates via other procedures. 
Our conclusions are presented in Section 7. 



2. Observations and target selection 

The data were taken using UVES (Ultra- Violet Echelle 
Spectrograph) on the European Southern Observatory Very 
Large Telescope (ESO VLT) at Cerro Paranal in Chile during 
22-24 July 2003, with cross-disperser #4. The spectra are in the 
wavelength range of 6100-9900 A, with a gap of ~80 A around 
8000 A due to the split between the two CCDs. Our observa- 
tions were taken at a time of superb seeing (~0.3" during the 
first 1 .5 nights, increasing to 0.6" for the remainder), allowing 
the use of the 0.3" slit. Thus, the spectra are of very high resolu- 
tion (R - 110, 000). Data were reduced using the UVES package 
for MIDAS and are shown in Fig.[T] 

Targets were selected from bright, red objects in the vicinity 
of galactic globular clusters and split into two categories based 
on the isocAM photometry of OFFR02: those with infra-red (IR) 
excess, {K - [12]o) > 0.65, referred to in this paper by the 'x' 
designator; and those without, {K - [12]o) < 0.65, referred to by 
the 'o' designator. The latter were selected from the 2MASS cat- 
alogue. A per-cluster summary of these objects can be found in 
Table[T] and a list of individual objects in Table|2] Radial veloc- 
ities taken from the spectra show all the stars to be likely cluster 
members, with the exception of the star observed in M54. The 



target selection is biased towards stars near the tip of the RGB 
and AGB, i.e. the most luminous red giant stars in the clusters, 
although we show that they do comprise an order of magnitude 
range in bolometric luminosity. Variability could have resulted 
in a bias towards stars in the bright maximum of their pulsa- 
tion cycle, but the IR amplitudes are only of order a tenth of a 
magnitude and well within the spread of the magnitudes. The 
inclusion of stars with known IR excess could in principle also 
bias the sample towards stars with higher mass-loss rates, but our 
analysis shows that stars without IR excess appear to experience 
similar mass-loss rates. 

In two cases in NGC362 - 2MASS J0103 1724-7050497 
(Bl) and 2MASSJ0103 1474-7050589 (B2) - the 2MASS data 
were not sufficient to determine which star was our target. In 
these cases, all possible counterparts (namely: o05a and o05b; 
and x02a, x02b and x02c, respectively) were also observed. 
Comparing the UVES acquisition and 2MASS images, we can 
see that the Bl is a blend comprising of two stars of similar 
magnitudes. We have assigned a flux based on Origlia et al.'s 
unpublished (J-K) data, split between the two stars on the basis 
of the comparative brightness in the acquisition images. The two 
stars have similar temperatures (see below), so we can expect 
them to have roughly similar (J-K) colours. B2 is a more com- 
plex blend of three bright stars (x02a, x02b and x02c) and two 
faint (unobserved) stars. The star x02c is more detached from 
the blend and has its own identifier (#1 1) in Valenti et al. (2004, 
hereafter VFO04): we do not consider it to contribute signifi- 
cant flux to the 2MASS blend. We thus assume that the major- 
ity of the flux observed in the 2MASS object is from the stars 
x02a and x02b and split the flux as before. These stars were also 
found to have similar temperatures. It thus appears likely that 
the 2MASS/ISOCAM identifiers are better associated with o05b 
and x02a, respectively. We show later that both sources appear 
to have higher mass-loss rates than their companions. 

Considerable literature already exists, both on the clusters 
and the individual stars we study. The cluster oj Centauri is well- 
noted for its chemically enriched sub-populations, with metal- 
licities spanning over an order of magnitude between [Fe/H] ~ 
-0.6 and [Fe/H] < -1.6 (Norris et al. 1996; Pancino et al. 2000; 
Sollima et al. 2005; van Loon et al. 2007) and helium abun- 
dances of up to 38% by fraction (Norris 2004). V42 is a well- 
known IR source (Origlia et al. 1995), studied recently spec- 
troscopically by van Loon et al. (2007) and with the Spitzer 
Space Telescope (Boyer et al. 2007, who find strong IR excess 
emission up to 24 pm). NGC 6388, one of the most metal-rich 
clusters, may also have a helium-enriched population (Moehler 
& Sweigart 2006). M 15 is a very metal-poor cluster but sign- 
ficantly the only cluster with a secure detection of ICM dust 
(Evans et al. 2003; Boyer et al. 2006) and gas (van Loon et al. 
2006a). Boyer et al. find M 15x01 is quiescent, but M 15x02 
to be mass-losing based on its IR colours. M54, a distant, 
but massive cluster with an unusually blue Horizontal Branch 
(Rosenberg et al. 2004) may be a component of the Sagittarius 
Dwarf Spheroidal galaxy (Ibata et al. 1994). NGC 362 is a small 
globular in a direction near that of the Small Magellanic Cloud 
(SMC). Finally, 47 Tuc, a more massive cluster near NGC 362 
and the SMC in the sky, has numerous variables which have re- 
cently been studied photometrically: Ramdani & Jorissen dis- 
cuss mass loss from V3 and VI 8, while Lebzelter & Wood 
(2005) discovered variability in several of our sample (LW8, 
LWIO, LW13 and LW19), and Ita et al. (2007) show V3 and 
LW13 to be on the giant branch tip, with V8 showing strong 
lljum excess. Spectroscopic observations in the mid-IR have 
been earned out by Lebzelter et al. (2006) on V3, V8 and VI 8, 
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Fig. 1. Overview of the VLT/UVES spectra of the target stars, sorted by spectral type (earliest at the top, latest at the bottom) and 
normalised to their respective 7000 A fluxes. The spectra have been smoothed for clarity over a 1000-pixel range (18-24A). Strong 
molecular absorption bands of TiO and VO are visible in the cooler stars and have been labelled, along with strong telluric features. 
Note the strong flux from M54x01 below ~7000 A when compared to the depth of its lines and molecular bands. 
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Table 1. List of clusters observed. 



NGC 


[Fe/H] 


Heliocentric 


Velocity 


Distance 


Absolute 


Stars with(out) IR Excess 


Ref- 


Number 




Radial Velocity 


Dispersion 




Visual 








erences 






(kms-') 


(km s-') 


(kpc) 


Magnitude 


With 


Without 


Unknown 




104 (47Tuc) 


-0.76 


-18.7 


11.5 


4.5 


-9.42 


9 








1 


362 


-1.16 


223.5 


7.5 


8.5 


-8.41 


3 


11 


2 


2 


5139 (coCen) 


-1.62t 


232.2 


19.6 


5.3 


-10.29 


1 








1 


6388 


-0.60 


81.2 


18.9 


10.0 


-9.42 


6 


12 





3 


6715 (M54) 


-1.58 


141.9 


14.2 


26.8 


-10.01 








1 


2 


7078 (M 15) 


-2.26 


-107.0 


11.0 


10.3 


-9.17 


2 








1 



Metallicity, radial velocity and distance values from Harris (1996). 

Other references: 1 - Paturel & Garnier (1992); 2 - McLaughlin & van der Marel (2005); 3 - Gebhardt et al. (1997). 
fMetal-rich ([Fe/H] — 0.6) and metal-intermediate (-1.3 < [Fe/H] < -1.0) populations are also believed to be 
present (Norris & da Costa 1995; Pancino et al. 2000; Sollima et al. 2005). 



and by van Loon et al. (2006b) on the above three stars and 
LWIO, with possible silicate emission present in V18, which 
may have recently undergone a thermal pulse (Lebzelter et al. 
2005). 



3. Stellar parameters 

3.1. Radial velocity 

The temperatures and radial velocities of our target stars were 
estimated using atlas9 model spectra (Kurucz 1993). A grid of 
model stellar spectra was created using a temperature range of 
3500-6000 K, in 250 K steps; metalHcity [Z/H] of -2.0 to -0.5, 
in steps of 0.5; and log(g) from to 1.5, in steps of 0.5 dex. 
In addition to the normal atomic line lists, we included titanium 
oxide bands (Kurucz 1999) to better replicate spectra at lower 
temperatures. 

A spectral type was assigned visually from standard refer- 
ence spectra (Pickles 1998), and an atlas9 template spectrum 
was chosen from the above grid that best replicated the tem- 
perature expected for that spectral type, the host cluster metal- 
licity, and an estimated gravity (based on Cioni et al. 2006a, 
2006b). Radial velocities were found by cross-correlating the at- 
LAs9 spectrum with the stellar spectrum in regions rich in atomic 
or molecular lines and clear of sky lines (avoiding the Ha and 
near-IR Ca II triplet lines). No significant difference between 
the velocities derived from the atomic- and molecular-line-rich 
regions was found. Spectra normalised in flux were also created 
using a running 1000-point (18-24 A) wide boxcar filter, chosen 
to minimise smoothing of molecular bands while maximising 
the removal of residual fringes that affect these red spectra. 

For the hotter objects, this yielded radial velocities accurate 
to within ~0.3 km s \ but for the cooler stars, where the strong 
molecular bands were not covered by the atlas9 models, the er- 
rors in radial velocity are much larger. In particular, those of 
47Tucx08, NGC 6388 ol2 and wCenxOl are over 10 km s"' 
(see Table |3j. Comparison to other data (van Loon et al. 2007) 
shows that the velocity for wCenxOl appears correct to within 
a few km s"'. This suggests that the formal errors we have cal- 
culated may be too large. 

The only object with a radial velocity not co-incident with 
its host cluster is that of M 54 xOl, suggesting it is either a field 
star or in a compact binary with a relatively massive compan- 
ion. The calculated luminosity for the distance of the cluster is 
not unreasonable for a red giant, and the spectrum appears to be 
that of a low gravity star. This suggests it may be in a binary 
system. However, attempts to cross-coiTelate the spectrum with 



5000 




K4 MO M2 
Spectral Type 

Fig. 2. Comparison of temperature determinations between two 
methods: open triangles represent temperatures derived from 
Kurucz's atlas 9 model fitting, filled circles represent 'cor- 
rected' temperatures (see text) derived from Fluks et al. (1994) 
and Levesque et al. (2006). The lines show the corrections ap- 
plied. The solid line represents Fluks et al.'s data for solar metal- 
licity, the dashed line represents the correction to [Fe/H] = -0.6 
(NGC 6388) and the dotted line to [Fe/H] = -0.76 (47Tuc). 
The data for NGC 362 at [Fe/H] = -1.16 have not been shown 
for clarity, as the correction is negligible in the observed range 
(bluewards of K5.5). The M5.5 outlier is wCenxOl, the fitting 
method for which is described in the text. 



the ATLAs9 models to find a binary component with a dififerent 
radial velocity were unsuccessful and its velocity suggests that 
it is most likely a field star, possibly in the Galactic Bulge (M 54 
is at I ^ 5.6°, -14.1°). 



3.2. Stellar temperature 

Stellar temperatures were calculated for the hotter objects via 
chi-squared minimisation using the atlas9 models, again using 
the host cluster metallicity and estimated gravity. This proved 
unsuitable for the cooler objects due to limitations of the atlas9 
models. For these objects (which typically have spectral class K4 
or cooler), temperature was estimated from their spectral type 
using the temperature-spectral type correlation from Fluks et al. 
(1994). These correspond to the 'uncorrected temperatures' in 
Tabled 
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Table 2. List of objects observed with positions, photometry and ancillary data. 



Our 


RA 


Dec 


J 


K 


IR 


Other 


Period 


Amplitude 


Notes 


ID 


(J2000.0) 


(J2000.0) 


Mag 


Mag 


Excess? 


Names 


(days) 


(K-band) 














NGC362 










oOl 


01 03 15.1 


-70 50 32 


9.93 


8.78 


N 


V16 


138 


0.5 


7 


o02 


01 03 13.6 


-70 50 59 


10.27 


9.51 


N 










o03 


01 03 13.8 


-70 51 09 


9.60 


8.79 


N 










o04 


01 03 12.6 


-70 50 38 


10.05 


9.13 


N 










o05a 


01 03 17.1 


-70 50 50 


9.91 : 


9.19 


N 








2,5 


o05b 


01 03 17.3 


-70 50 50 


9.95 : 


9.27 


: N 








2,5 


o06 


01 03 13.6 


-70 50 37 


9.68 


8.74 


N 






0.15 : 


7 


o07 


01 03 20.1 


-70 50 55 


9.58 


8.75 


N 










o08 


01 03 10.7 


-70 50 54 


9.85 


8.89 


N 










o09 


01 03 10.9 


-70 50 59 


10.66 


9.83 


N 










olO 


01 03 14.7 


-70 51 15 


10.38 


9.67 


N 










xOl 


01 03 19.0 


-70 50 52 


9.96 


8.96 


Y 










x02a 


01 03 14.7 


-70 50 59 


9.64: 


8.27 


: Y? 








2 


x02b 


01 03 14.9 


-70 51 00 


9.82: 


8.45 


: Y? 








2,5 


x02c 


01 03 14.5 


-70 50 59 


10.45 


9.64 


Y 








3,5 


x03 


01 03 13.7 


-70 51 14 


10.10 


8.82 


Y 




















47Tuc 










xOl 


00 24 02.6 


-72 05 07 


7.73 


6.75 


Y 


LWIO 


110 






x02 


00 24 07.8 


-72 05 09 


8.80 


7.79 


Y 










x03 


00 24 09.4 


-72 04 49 


8.37 


7.35 


Y 






0.15 : 


7 


x04 


00 23 57.7 


-72 05 30 


8.77 


7.75 


Y 


LW8 


27 


0.15 : 


7 


x05 


00 24 23.2 


-72 04 23 


8.30 


7.16 


Y 


LW19 


40 






x06 


00 24 08.6 


-72 03 55 


8.09 


6.84 


Y 


V8? 


155 


0.3 


7 


x07 


00 24 07.9 


-72 04 32 


7.82 


6.56 


Y 


LW13 


65 






x08 


00 25 16.0 


-72 03 55 


7.66 


6.45 


Y 


V3 


192 


0.6 


1,7,8 


x09 


00 25 09.2 


-72 02 40 


8.67 


6.56 


Y 


V18 


83 




1 












NGC6388 










oOl 


17 36 14.9 


-44 43 12 


10.73 


9.31 


N 


V9 








o02 


17 36 11.4 


-44 44 08 


10.02 


8.80 


N 


V12 




0.25 : 


7 


o03 


17 36 20.2 


-44 44 32 


11.01 


9.65 


N 










o04 


17 36 17.6 


-44 43 46 


10.49 


9.29 


N 










o05 


17 36 16.7 


-44 43 28 


10.86 


9.64 


N 










o06 


17 36 17.7 


-44U12 


10.37 


9.27 


N 








1 


o07 


17 36 15.5 


-44 44 26 


10.72 


9.31 


N 










o08 


17 36 15.3 


^4 44 31 


10.91 


9.47 


N 










o09 


17 36 18.9 


-44 44 14 


11.31 


10.13 


N 










olO 


17 36 16.0 


-44 43 49 


10.73 


9.72 


N 










oil 


17 36 20.4 


-44 44 16 


11.12 


9.97 


N 










ol2 


17 36 18.9 


-44 44 04 


11.09 


9.81 


N 










xOl 


17 36 18.1 


^43 25 


10.80 


9.75 


Y 


V8? 




0.4 


7 


x02 


17 36 13.1 


^43 05 


10.65 


9.47 


Y 










x03 


17 36 15.1 


^43 33 


10.27 


8.99 


Y 


V3 


156 


0.3 


7 


x04 


17 36 21.4 


^43 42 


11.15 


9.87 


Y 










x05 


17 36 18.8 


-44 43 46 


10.60 


9.25 


Y 










x06 


17 36 16.0 


-44 44 36 


10.64 


9.19 


Y 




















M15 










xOl 


21 29 58.9 


+12 10 18 


11.18 


10.50 


Y 


AC761 








x02 


21 29 58.2 


+12 09 47 


10.63 


9.92 


Y 


AC414 






9 












M54 










xOl 


18 55 04.5 


-30 29 35 


11.74 


11.12 


Y 




















wCen 










xOl 


13 26 46.4 


-41 29 30 


9.70 


8.75 


Y 


V42 


149 


0.5 


7,8 



Notes: Astrometry/photometry from OFFR02 or their unpublished data unless noted otherwise. 

Photometry sources: (1) 2MASS (corrected to the photometric system of OFFR02); (2) UVES acquisition images 



(see text); (3) VFO04. Colons denote uncertain values. 
Astrometry sources: (4) 2MASS; (5) UVES acquisition images (see text); (6) VFO04 

Variability sources: LW indicates variables from Lebzelter & Wood (2005). (7) Matsunaga et al. (2007), Ita et al. (2007) 

or Matsunaga (private communication); (8) Feast et al. (2002). 
(9) identified with 2MASS J21295815+1209466, which is bright at 24jum (M24 =1.66 - Boyer et al. 2006). 
AC identifiers are for (non-variable) sources from Auriere & Cordoni (1981). 



6 



I. McDonald & J. Th. van Loon: Dust, pulsation, chromospheres and mass loss in red giants in Galactic GCs 



4600 



4400 - 



4200 



ra 4000 



o 
o 
O 



3800 
3600 



3400 



3200 



3000 



2800 



~i 1 1 1 1 1 1 r 



2800 3000 3200 3400 3600 3800 4000 4200 4400 4600 
Corrected Temperatures (K) 

Fig. 3. 'Corrected' temperatures from spectral typing or at- 
LAs9 model fitting versus temperatures determined from (J-K) 
colour using Houdashelt et al. (2000). The two methods appear 
poorly correlated, as discussed in the text. Key: NGC 362: plus- 
signs, 47 Tuc: crosses, NGC 6388: asterisks, M 15: open squares, 
oj Cen: filled square. 



The Fluks et al. relations are valid for solar metallicity stars. 
As metallicity affects the value of eff'ective temperature deter- 
mined for a given spectral type, it was necessary to correct the 
cool objects' temperatures for the substantially sub-solar metal- 
licities of the clusters. To do this, we used a similar relation for 
LMC, SMC and Galactic red supergiants from Levesque et al. 
(2006) and extrapolated the metallicity dependence at each spec- 
tral sub-type out to our clusters' metallicities using a second- 
order polynomial interpolation. We assumed [Z/H]lmc = -0.3 
and [Z/H]sMC = -0.6 (Russell & Dopita 1990), and [Z/H]mwg = 
0, and used a linear interpolation between sub-types where nec- 
essary. Overall, this has the effect of lowering the temperature of 
the stars by around 100-200 K. 

This approach appears to be satisfactory in all cases except 
the M-giant in u) Cen, for which this method would yield a tem- 
perature of 2000 K. Here, a linear extrapolation using the Milky 
Way and SMC data at spectral type M2 (the last datapoint in the 
Levesque et al. data) was used instead of a parabola. Due to the 
uncertainty of the extrapolation and possible errors introduced 
by this correction process, we suggest that the temperature we 
find is still rather low. Were this star part of the metal-rich sub- 
population, its temperature would be significantly higher (up to 
3250 K, at [Z/H] = -0.6). We have, however, no evidence to 
support this theory, and van Loon et al. (2007) find that the M- 
type stars in o) Cen may have metallicities similar to the cluster 
average. 

As the Levesque et al. data are for supergiants and not for 
the less massive and less luminous giants, we have to further 
correct for the difference. We perform a crude correction, using 
the offset between the temperatures of Levesque et al.'s super- 
giant galactic sample and Fluks et al.'s red giant (galactic) data. 
For each particular subtype, we subtract the difference from our 
result; again we interpolate where necessary. This has an effect 
of lowering the effective temperatures of types M0-M3 by ~ 100 
K, but little difference outside this range. These values corre- 
spond to the 'corrected temperatures' in Table [3] They are plot- 
ted alongside the atlas9 model fits in Fig.|2] 

A comparison between the two methods for deriving temper- 
ature shows that the atlas9 models are still consistently around 



200 K cooler than the corrected temperatures from spectral typ- 
ing (Fig.|2|, suggesting our corrections have not been as large as 
necessary. This may be due to the uncertainty in extrapolating 
the empirical corrections to low metallicities. One would pre- 
sume this offset still exists at cooler temperatures where the at- 
LAs9 models have not made a fit, however as this is uncertain 
we have not made any further corrections for this in our results. 
This suggests that metal-poor spectra may correspond to even 
lower temperatures (around the lower error values) than calcu- 
lated here. 

Comparisons were made with the colour-temperature tables 
in Houdashelt et al. (2000), using (J-K) colour to estimate a tem- 
perature. The temperatures derived using this method are plotted 
against the temperatures found using the above methods (after 
corrections) in Fig. [3] The Houdashelt et al. relations are cali- 
brated only for temperatures > 4000 K, thus it is not surprising 
that we observe a poor correlation between the two methods. 

3.3. Physical Parameters 

Physical parameters of the stars have also been calculated (Table 
[3]l. Luminosities were found from K-band magnitudes and us- 
ing bolometric corrections made from Montegriffo et al. (1998) 
using a linear interpolation for metallicity. Radii were then com- 
puted using the 'corrected temperatures' described above and, 
using these, escape velocities were also computed. A mass of 
0.8 Mq was assumed for all stars, though some stars near the 
AGB tip may have lost up to ~30% of this, which would lower 
their escape velocities by up to ~17%. 

3.4. Analysis 

We display a Hertzsprung-Russell diagram in Fig. ID includ- 
ing isochrones at 12.6 Gyr from Cioni et al. (2006a, 2006b). 
There is significant scatter in the plots, likely due to the uncer- 
tainty in temperature of the stars, though the stars of metal-rich 
and metal-poor clusters clearly match the metal-rich and metal- 
poor isochrones, respectively. The position of the RGB tip, at 
around 2000 Lq, does not change appreciably with metallicity 
nor within the likely range of ages. Assuming, therefore, that 
the distances to the clusters are correct, we expect that the stars 
above the RGB tip are on the AGB (below the tip it is difficult to 
tell which branch they belong to) suggesting more than half of 
the sample is on the AGB. Note that, due to their intrinsic vari- 
ability, the stars will naturally display a range in both luminosity 
and temperature: a 0.5 (bolometric) magnitude variability would 
change the luminosity by a factor of 1.58, and the temperature 
by about 130 K (assuming the net flux at the surface is roughly 
constant). 

A period-magnitude diagram is shown in Fig.|5] using peri- 
odicity data from the Updated Third Catalogue of Variable Stars 
in Globular Clusters (Clement 1997) and from Matsunaga et al. 
(2007), Ita et al. (2007) and Matsunaga (private communication). 
Frequency modes from Ita et al. (2004) have also been shown. 
These relations are for the LMC, so we expect our more metal- 
poor sample to be slightly brighter for a given period (Wood 
1990). It is notable that very few stars in our sample have known 
variability data, and we discuss this in Sect. 6.1. 



4. Ha and Ca II profiles 

In an effort to find mass-losing stars, we investigated the Ha 
(6562.8 A) and calcium IR triplet lines (8498, 8542 and 8662 
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ID 


Spectral * 


Uncorrected 


Corrected 


Colour-Derived 


Radial 


Luminosity 


Radius 


Escape 




Type 


Temperature 


Temperature 


Temperature 


Velocity 






Velocity 






(K) 


(K) 


(K) 


(km s-') 


(Lo) 


(R0) 


(km s-') 



oOl 
0O2 
o03 
o04 
o05a Kl 
o05b Kl 
0O6 



o07 
0O8 
o09 
olO 
xOl 



K5.5 
KO 
Kl 
Kl 



K5.5 
K3 
K2 
KO 
KO 
Kl 



x02a K3 

x02b K3 

x02c KO 

x03 K2 



NGC362 



F 
K 
K 
K 
K 
K 
F 
K 
K 
K 
K 
K 
K 
K 
K 
K 



3954 +'7 

4200 

3900+1^ 

4125 

4100 

tiw 250 

4125 +2^" 
-230 

3Q^4 +' " 

3950 :||| 

3850 _25o 
4350+250 

4600+?!^ 



4000 
3950 
3950 
4375 
3950 



3867!]^, 
4200 +|i! 
3900 +|g 
4125 *™ 

250 

4100 +250 

tiw 250 

4125 
3867!!^. 
3950 :||| 

3850 _25o 
4350+250 
4600 +2| 
4000 tg^ 
3950+1° 
3950 +2| 
4375 +2| 
3950+2° 



3337 


2.5 ± 


4.3 


2415 ± 


102 


109 ± 


5 


53 ±6 


4191 


9.3 ± 


1.0 


988 ± 


86 


59 ± 


6 


72 ± 10 


4055 


-16.7 ± 


1.0 


2327 ± 


226 


105 ± 


12 


54±8 


3782 


-A.2± 


0.3 


1468 ± 


131 


75 ± 


8 


64±9 


4527 


-6.8 ± 


0.3 


2103 ± 


189 


91 ± 


10 


58 + 9 


4527 


-5.8 ± 


0.3 


2392 ± 


213 


95 ± 


10 


57 ±8 


3732 


0.7 ± 


2.7 


2505 ± 


105 


111± 


6 


52 ± 6 


4003 


-8.6 ± 


0.3 


2327 ± 


222 


103 ± 


12 


54±8 


3682 


7.3 ± 


0.3 


2182 ± 


215 


105 ± 


12 


54±8 


4003 


10.3 ± 


0.3 


665 ± 


54 


45 + 


5 


82 ± 12 


4343 


1.7 ± 


0.4 


659 ± 


48 


40 ± 


4 


87 ± 12 


3589 


-1.6 ± 


0.3 


1866 ± 


175 


90 ± 


10 


58 ±9 


3232 


-1.4 ± 


0.3 


2599 ± 


248 


109 ± 


12 


53 ±8 


3232 


6.2 ± 


0.4 


2824 ± 


269 


113± 


13 


52± 8 


4239 


0.9 ± 


0.3 


2027 ± 


164 


78 ± 


8 


62 ± 9 


3220 


^.1± 


0.3 


2182 ± 


208 


99 ± 


11 


55 ± 8 













47Tuc 


















xOl 
x02 


M4 
K4 


F 
F 


3574 +'2 
4147 +f 


3351 +'2 
4172+^^ 


3724 
3659 


-7.4 ± 
7.2 ± 


6.4 

2.4 


6172 ± 
1375 ± 


427 
51 


232 + 
71 + 


19 

3 


36 + 
66 + 


5 
7 


x03 


K5 


F 


3993 +" 


3985 +" 


3639 


-10.5 ± 


0.6 


2325 ± 


103 


101 + 


5 


55 + 


6 


x04 


K4 


F 


4147 +1 


4172+1 


3639 


-5.6 ± 


3.1 


1427 ± 


53 


72 + 


3 


65 + 


7 


x05 


MO 


F 


3895 
3371 +15^ 


3777 *f 
3144+1^2 


3427 


10.4 ± 


3.0 


3180± 


135 


131 + 


6 


48 + 


5 


x06 


M5.5 


F 


3299 


-2.7 ± 


7.3 


6463 ± 


448 


270 + 


22 


34 + 


4 


x07 


M4 


F 


3574 +92 


3351 +''2 

jjji -140 


3290 


4.5 ± 


6.4 


7352 ± 


509 


254 + 


21 


35 + 


4 


x08 


M7 


F 


3188 +{25 
4147 +^2 


2982 +'25 

-125 

4172 +^2 


3340 


-6.5 ± 


11.7 


10338 + 


711 


380 + 


31 


28 + 


4 


x09 


K4 


F 


2988 


10.8 ± 


2.6 


4270 + 


158 


125 + 


5 


49 + 


5 













NGC 6388 
















oOl 


M4 


F 


3574 +'2 


3545 


3578 


4.4+4.1 


1820 + 


117 


121 + 


9 


50 + 


6 


0O2 


M3 


F 


3666 +™ 
4050 +2'o 

_25o 


4050+1 


3328 


-16.0+4.0 


4055 + 


171 


168 + 


8 


43 + 


5 


o03 


K3 


K 


3218 


23.8 + 0.3 


1331 + 


122 


74 + 


8 


64 + 


10 


o04 


MO 


F 


3895 +»» 
3993 *ll 


3808 +1 
3999 *ll 


3352 


31.7+3.0 


2168 + 


79 


107 + 


5 


53 + 


6 


o05 


K5 


F 


3328 


3.4+3.1 


1380 + 


54 


77 + 


3 


63 + 


7 


0O6 


Ml 


F 


3810 


3735 t'jl 


3488 


3.1 + 2.5 


2312 + 


77 


115 + 


4 


52 + 


6 


o07 


K4 


K 


4400 +250 
-250 


4400 +250 

_250 


3194 


7.2+0.8 


1445 + 


115 


65 + 


6 


68 + 


10 


0O8 


M2 


F 


3736 +^* 


3663 


3185 


-19.5 + 3.8 


2032 + 


66 


112 + 


4 


52 + 


6 


o09 


M1.5 


F 


3773 +™ 


3699 !™ 


3376 


5.9+3.0 


1077 + 


35 


80 + 


3 


62 + 


7 


olO 


K5.5 


F 


3954!]' 


3954 +■'■' 

jyj^ -101 


3659 


-16.4 + 2.6 


1331 + 


54 


78 + 


4 


63 + 


7 


oil 
0I2 
xOl 


M3.5 

M4 

M3.5 


F 
F 
F 


3620! , 
3574 +'- 

-'-^'^-140 

3701 +^j 


3513! , 
3428+1 
3585 gj 


3411 
3273 
3578 


0.1 + 8.4 
-12.5+ 11.3 
-0.1 + 4.4 


1406 + 
1722 + 
1644 + 


74 

111 

61 


101 + 
117 + 
105 + 


6 
9 
5 


55 + 
51 + 
54 + 


6 
6 
6 


x02 


Ml 


F 


3810 +^J 


3735 


3376 


-5.7+2.4 


1923 + 


64 


104 + 


4 


54 + 


6 


x03 


M1.5 


F 


3773 +™ 
3666 +™ 
3773 +™ 
3666 !™ 


3699 *ll 


3273 


-19.6+3.4 


3076 + 


101 


135 + 


5 


48 + 


5 


x04 


M3 


F 


3545 +™ 


3273 


26.4 + 4.0 


1514 + 


64 


103 + 


5 


54 + 


6 


x05 


M1.5 


F 


3699 !| 


3224 


7.9+2.4 


2421 + 


79 


119 + 


5 


51 + 


5 


x06 


M3 


F 


3545 


3182 


9.2+3.6 


2831 + 


119 


141 + 


7 


47 + 


5 












M15 
















xOl 

x02 


G8 
G8 


K 
K 


4500 !25« 
4400 !25° 


4500 !25° 
4400 !25° 


4448 
4343 


-14.7 + 7.7 
10.3 + 5.1 


489 + 
891 + 


37 
71 


36 + 
51 + 


3 
5 


92 + 
77 + 


13 
11 












M54 
















xOl 


K2 


K 


3800 !25o 


3800 !25o 




-182.9 + 0.2 
























a>Cen 
















xOl 


M5.5 


F 


3371 !;32 


2675 !^^2 


3707 


44.3 + 26.2 


2111 + 


161 


213 + 


19 


38 + 


5 


Notes: Column 3: F 


- from spectral typing, using 


; Fluks et al. (1994); K - derived from Kurucz's atlas9 models (see 



Section 3). The colour-derived temperature is from (J-K) colours, from Houdashelt et al. (2000) (see text). Radial velocity 
is relative to the mean motion for the cluster from Harris (1996), listed in Table[T] See text for details on M54x01 and 
aiCenxOl. Formal errors have been assigned on the basis of a 0.1 dex uncertainty in metallicity, a 100 pc uncertainty in 
distance and a 0.1 M^ error in mass. 
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Fig. 4. Hertzsprung-Russell diagrams of our targets. Top panel: 
stars identified by cluster: NGC 362 - plus signs, 47 Tuc - 
crosses, NGC 6388 - asterisks, M 15 - boxes, M54 - filled 
box, a>Cen - circle. Bottom panel: IR excessive stars (crosses) 
and non-excessive stars (boxes). The temperature for our w Cen 
target has probably been underestimated. Also shown are 
isochrones for [Z/H] - -0.67, -1.28 and -1.68 (labelled on top 
panel) from the models of Cioni et al. (2006a, 2006b) at an age 
of 12.6 Gyr. 



A) for evidence of asymmetry and velocity shift. These lines are 
shown in Fig. |6] To accomplish this, we calculated line bisec- 
tors by interpolating the line profiles to uniform steps in flux. At 
each given flux, the two points at either side of the line centre 
(as determined by the line's rest wavelength with respect to the 
calculated radial velocity) were averaged. 

Examination of the calcium triplet line bisectors shows most 
to be fairly symmetrical. Some show slight blue- or red-shifted 
cores of the order 1-3 km s"' , but these are not usually visible in 
all three lines and can probably be attributed to statistical noise, 
rather than as certain indicators of mass loss or pulsation (the 
exceptions being 47Tucx04 and x09, which show 4-5 km s ' 
blue-shifted cores in all three lines). We visually estimate that 
the uncertainty on the shifts is of the order of 0.5-2 km s de- 
pending on the signal-to-noise of the data and eff'ects of intrud- 
ing metal lines. This is similar to, or less than, the radial velocity 
uncertainties. The Ha line shows similar features in some cases, 
but overall shows much more prominent emission and more sig- 
nificantly blue-shifted absorption at velocities of 2-12 km s"'. 
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Fig. 5. Period-magnitude diagram for stars with variability data. 
The lines show the fundamental mode (longest period), and first 
and second overtone (harmonic) periods, using the relations for 
the LMC given by Ita et al. (2004). The corresponding Ha pro- 
files are plotted near the positions of the stars. Clearly, Ha emis- 
sion is strongest in the fundamental mode and is therefore most 
likely linked to shocks in the pulsating atmosphere. 



The Ha lines also tend to be very deep, with over 90% absorp- 
tion in the core in some cases. 

In both elements, the lines were predominantly found to be 
'banana-shaped', with a reddened kink at approximately 30- 
50% of the peak absorption depth. As such, the profiles are not 
well suited to the usual method of finding a 'velocity shift'; thus, 
we calculated the radial velocities of these 'kink' points and of 
the line cores (negating anything appearing to be statistical scat- 
ter due to Poisson noise) with respect to the stellar velocities 
described above. The results are included in TablelH 

Notable exceptions to the above features appear in gen- 
erally more luminous and cooler objects. Strong emission is 
clearly present in 47Tucx08 and wCenxOl in Ha and (to a 
lesser extent) in Ca II. Lesser versions of this can also be seen 
in NGC362o01 and to some extent in NGC 6388x03 (Fig. 
5). These four stars show a much larger K-band variability 
(hence much stronger pulsations) than the remainder of the ob- 
jects for which we have variability data (see Table |2]). This 
emission is therefore almost certainly due to heating by shock 
waves caused by pulsations. NGC362o01 and wCenxOl have 
also been shown to have variable Balmer line emission (Lloyd 
Evans 1983a, 1983b), further suggesting this is pulsation related. 
47Tucx09 also shows variable Ha absorption (Lloyd Evans 
1984). 

Furthermore, it is interesting to note that the two most ex- 
treme objects show qualitatively very similar profiles in Ha and 
to a lesser extent in the calcium triplet. This is particularly inter- 
esting given their clearly different metallicities: 47 Tuc x08 is es- 
timated to be at 3-8 times of the metallicity of cj Cen xO 1 , which 
is visibly manifest in the other, stronger spectral lines of the for- 
mer in Fig.|6l This has obvious implications for the variance (or 
otherwise) of mass-loss rate with metallicity. 

Weaker emission is also apparent in warmer stars. This has 
been noted for many giant stars in the literature and two re- 
lated explanations have been put forward to model this. Cohen 
(1976) assumes a cooling circumstellar envelope surrounding 
the star at around two stellar radii. Another theory, by Dupree 
et al. (1984), states that Ha emission could arise from chro- 
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k [A] \ [A] A [A] X [A] 

Fig. 6. Ha and near-IR calcium triplet line profiles of the target stars, sorted by spectral type as in Fig. [T] and normalised to their 
continuum fluxes at the wavelengths of the respective lines. The spectra have also been shifted to laboratory wavelengths. Note 
the variety of asymmetric emission and absorption in the lines, especially in Ha. We consider the strong Ha emission from the 
two latest-type objects to be shock emission from pulsation. The similarity of the profiles is remarkable, despite very different 
metallicities. 
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mospheric emission from a similar, hot shell source. Cacciari 
et al. (2004) have suggested that active chromospheres are the 
dominant producer of line emission in these stars. Mauas et al. 
(2006) then used these chromospheric lines to model RGB mass 
loss in good agreement with evolutionary theory and previous 
estimates (Dupree et al. 1990, 1994) - notably they find the need 
for an outward velocity field in the chromosphere, which we also 
find (see discussion). 

It is an obvious suggestion, backed up by more detailed mod- 
elling by Dupree et al. and Mauas et al., that blue-shifted absorp- 
tion cores are evidence of outflow from the stars. The fact that 
we see red-shifted cores in only two cases (two of the pulsation- 
shocked objects) suggests that there is little or no inflow to the 
stars, or that this only occurs at very particular phases of the pul- 
sation cycle. We assume that there is therefore a net mass loss 
from the sample. 

5. A simple model for estimating mass loss 

5. 1 . The model 
5.1 .1 . General form 

In our attempt to quantify mass-loss rates from our sample of 
stars, we have created a simplistic model, based on some crude 
approximations, which we use to model the Ha line profile. 
Estimations of mass-loss rates based solely on line profiles are 
known to be notoriously uncertain, in part due to their inherent 
degeneracy in parameters such as velocity and turbulence (see 
e.g. Hagen 1978; Hagen et al. 1983). We address these problems 
in Section |572l 

As illustrated in Fig. |7] we assume a spherically-symmetric 
wind expanding with constant speed Vm, emanating from a non- 
rotating star. We assume that the emission from the surface of 
the star in the region of the Ha line can be simplified to a flat 
continuum, modified by a Gaussian-shaped absorption line, cen- 
tred on Ha. This line profile is then modified by the stellar wind 
and finally convolved with a Kurucz stellar model with the same 
metallicity, temperature and gravity as the star, broadened by a 
macroturbulence parameter Vtmb, but without any hydrogen lines. 
Calculating the model in this way, rather than starting with the 
ATLAs9 model as the output from the photosphere allows an ana- 
lytical representation of the photospheric emission which greatly 
improves calculation speed. 

The width of the photospheric Ha line is given by: 



w 



— + 



tuib 



(r = R,), 



(1) 



where k is Boltzmann's constant, Tph is the effective temperature 
of the photosphere, jUh is the (reduced) atomic weight of hydro- 
gen and v'rjns is the root mean square velocity due to thermal 
motions (additional contributions to motion in the photosphere 
are covered by the 'turbulence' term, vturb)- Thus, the flux sub- 
tracted from the continuum by the photosphere, in terms of the 
continuum flux, can be described as simply; 



(2) 



where v is the radial velocity relative to Ha and Tph is the optical 
depth of the photospheric absorption in the line centre. 

When calculating the wind contribution, we only consider 
the absorption and emission from bound-bound transitions for 
simplicity - unless there is a chromosphere or other temperature 




Fig. 7. Illustration of our wind model, defining the symbols used 
in the text. 



enhancement (discussed later in this section), the temperatures 
in cool giant winds are never high enough to create sufficient 
amounts of free electrons that other contributions to opacity and 
line broadening become important. Thus, the observed flux from 
the star in velocity interval dv can be modelled as: 



F(v)dv _ ^ Apbiv)dv ^ E^ind{v)dv 



i{v)dv 



Fctmdv 



Fctmdv 



Frimdv 



Fctmdv 



(3) 



where Fcm is the continuum emission near Ha; E„\nd is the emis- 
sion component of the wind, due to photons emitted by sponta- 
neous emission directed into the line of sight; and A„ind is the ab- 
sorption component of the wind, due to photons absorbed/rom 
the line of sight radiation from the stellar surface. 

5.1.2. The absorption term 

The absorption from the column in front of the star cannot be 
adequately modelled by a simple column of absorbing material, 
but must take into account the spherical nature of the star. Thus, 
we define an impact parameter, b, which represents the physical 
distance in the plane of the sky from the centre of the star. 

We define the number density of absorbers (i.e. hydrogen 
atoms in the n = 2 level for Balmer a) at a given point in the 
wind as being n2ir). This can be found in local thermodynamic 
equilibrium (LTE) through the Saha-Boltzmann equations (the 
validity of an LTE environment is discussed later). Using this, 
we can calculate the optical depth integrated over all velocities 
in a particular radius bin dr as being dT{r) = n2{r)o-dr for the 
velocity-integrated absorption cross-section cr. Here, we use: 



ne 
mc 



(4) 



giving a value of 6.537 X 10^'^ m^ using Gray (1992). 

Integrating in annuli (b) over the observed stellar surface and 
outwards in radii, we can define the integrated optical depth of 
the absorption column at velocity v as: 



T(V) 



1 



y/TjTwir) 



exp 



-(v + v^Rjby 

2wHr) 



Inbdb 
X —dT(r)dr, 



(5) 
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Table 4. Radial velocities of the Ha and Ca II lines as described in Sect. 4. 



ID 


Hot 


Ha velocity 


Ca II velocity 


Ca II velocity 


Ca II velocity 




profile 


6563 A 


8498 A 


8542 A 


8662 A 




appearance 


Core 


Kink 


Core 


Kink 


Core 


Kink 


Core 


Kink 








NGC 362 












oOl 


MOL, BE, RE 


3.3'' 


-3.1'" 


4.0 


10.5 


-0.9 


10.9 


-0.6 


6.7 


o02 


PHOT, be 


-1.9 : 


4.8" 


2.7'' 


0.6'' 


2.2* 


2.2* 


1.9* 


1.9* 


o03 


PHOT, be 


0.7 


5.0" 


1.7 


3.5 


1.3 


4.0 


1.3 


3.9 


o04 


PHOT, be 


-2.4 : 


6.5° 


0.5 


1.7 


0.7 


2.7 


0.1 


2.1 


o05a 


PHOT, re? 


-6.5 


-0.4 


1.5* 


1.5* 


0.9* 


0.9* 


0.0 


2.0 


o05b 


PHOT, wide 


-1.8 : 


0.9" 


2.0* 


2.0* 


0.0 


2.0 


0.5 


2.9 


o06 


phot, BE, RE 


-1.6 : 


2.5" 


2.8* 


2.8* 


1.0 


2.2 


1.3 


4.2 


o07 


PHOT, be, re 


-6.2 


4.1" 


3.5* 


3.5* 


-0.8 


4.8 


3.3* 


3.3* 


o08 


PHOT, re 


-7.9 : 


-0.1 


1.3 


1.8 


0.8 


2.0 


1.5 


2.5 


o09 


PHOT, be, re 


-5.4 : 


2.0" 


-0.5 


1.4 


-1.9 


2.6 


1.0 


1.5 


olO 


PHOT 


2.9 


9.3" 


2.1 


5.1 


2.9 


6.0 


2.5 


6.4 


xOl 


PHOT, be?, re 


-3.8 : 


0.4" 


2.0* 


2.0* 


-2.0 


2.9 


0.4 


4.1 


x02a 


PHOT, asy 


-10.1 


2.6° 


-1.8 


1.5 


-5.4 


1.3 


-2.1 


2.2 


x02b 


PHOT, asy, re 


-12.5 


-1.7 


1.0* 


1.0* 


-4.3 


-0.0 


-1.5 


0.9 


x02c 


PHOT 


-1.5 : 


-3.3" 


1.5* 


1.5* 


0.9* 


0.9* 


0.6 


2.0 


x03 


PHOT, RE 


-3.5 


1.0 


1.4'' 


-2.1'' 


-7.2 


-0.1 


0.0 


-2.1 








47Tuc 












xOl 


MOL 






0.7 : 


2.8 


0.8 :* 


0.8 :* 


2.0 


0.0 




x02 


PHOT, asy 


-8.8 


-0.9 


-0.9 : 


1.2 : 


-6.7 


-0.1° 


-1.0 


0.5 




x03 


PHOT, narrow 


-1.0 


4.4 


4.4'' 


0.4'' 


-1.9 


4.2 


0.0 


3.0 




x04 


PHOT, asy 


-6.5 : 


-1.9° 


-4.1 


1.5 


-6.2 


2.7 


-7.0 


0.5 




x05 


PHOT, mol 


-4.8 


4.0° 


-1.0 


1.7 


-2.6 


1.8 


-1.5 


1.0 




x06 


Noisy continuum 










-0.5 : 


1.0 :° 


0.0 


1.0 


a 


x07 


Noisy continuum 






6.0 :'' 


0.1 :'' 


4.7* 


4.7* 


-1.0 


1.0 


a 


x08 


BE, RE 


-5.0 : 


-2.0 :" 


1.0 : 


3.0 


0.0 


1.5 


0.3 


1.0 


a 


x09 


PHOT, asy 


-7.8 


-0.9° 


-4.6 


-0.3 


-7.9 


-1.2 


-8.0 


-1.5 


a 








NGC 6388 












oOl 


PHOT, mol 


-6.7 


0.6 


-0.4 


1.6 


-3.9 


1.6° 


-1.3 


1.4 


o02 


phot 


4.7 


9.5 


8.5'' 


-0.4 :'' 


3.1 


7.4 


7.0 


1.4 


d 


o03 


PHOT 


-4.4 


0.9 


-1.7 


1.2 


-0.6* 


-0.6* 


-2.7 


0.5 


o04 


phot 


-2.3 : 


0.1 


0.1* 


0.1* 


-2.8 


-0.6" 


-1.6 


-0.1 


o05 


PHOT, asy 


-9.7 


2.6" 


1.0* 


1.0* 


-2.3 


0.6 


-1.8 


1.0 


o06 


PHOT, mol 


-6.7 


-0.5 


-0.1* 


-0.1* 


-5.0 


-0.3" 


-3.3 


-0.7 


o07 


phot, mol, asy 


-5.7 


1.9 


-0.6 


1.0 


-3.5 


-0.9 


-2.3 


-0.3 : 


o08 


phot, mol 


-2.3 


0.6 


-0.5 


1.8 


-1.0 


0.8 


-0.2 


1.2 


o09 


PHOT, wide 


-0.0 : 


2.0 


1.7'' 


-0.1'' 


1.1* 


1.1* 


1.2* 


1.3* 


olO 


PHOT, asy 


-7.4 


-2.6° 


4.5'' 


2.1 :'' 


-0.2 


3.3 


1.9* 


1.9* 


oil 


PHOT, be? 


1.0* 


1.0* 


0.7* 


0.7* 


-0.1* 


-0.1* 


0.8* 


0.8* 


ol2 


PHOT, mol 


-6.3 


-0.6 


0.0 


1.6 


-2.6 : 


-0.8" 


-2.5 


0.7 


xOl 


phot, BE, RE 


-0.9 


1.9 


4.1'' 


0.4 :'' 


3.7'' 


-3.4'' 


3.3 


-1.6'' 


x02 


PHOT 


-3.7 


2.6 


0.2* 


0.2* 


-1.8 


0.9 


-0.8 


1.2 


x03 


BE 


24.3" 


19.4" 


17.0'' 


1.9 :'' 


11.7'' 


-4.7 :'' 


14.7'' 


-0.2 


d 


x04 


phot, mol 


-2.3 : 


2.9 


0.4* 


0.4* 


0.5* 


0.5* 


0.6* 


0.6* 


x05 


PHOT, mol, asy 


-3.3 


6.0° 


-0.5 


2.8 


-3.3 


1.7 


-1.2 


3.3 


x06 


MOL 


1.0 :* 


1.0 :* 


2.8 :'' 


-0.4 :'' 


-2.6 


0.3 


0.5* 


0.5* 








M 15 












xOl 


PHOT 


-4.2 


3.3" 


0.1'* 


-1.7'' 


0.0'' 


-3.6'' 


-0.1* 


0.1* 


x02 


PHOT, be, re? 


-12.9 


-0.5° 


1.2 


2.7 


0.8 


2.4 


-1.7 


3.9 








wCen 












xOl 


BE, re 


-3.0" 


-4.0" 


+1.0" 


+0.5" 


+0.5" 


-1.8" 


0.0" 


-4.0" 



Positive values denote red-shifts (in km s '). Colons denote ill-defined values due to noise. Abbreviations: 
phot - photospheric profile visible; asy - asymmetric; mol - molecular bands; be/re - blue/red emission 



wing; capitalisation indicates strong features. Notes: " kink not present, value given near continuum level; 
* no significant velocity shift in bisector, approximate average given; ' core at a greater red-shift than kink 
or no kinlc, due lo an emission line or ('') unusual prolile. 
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with VooR*/b being the projected velocity of the wind into our 
line of sight. We can thus define the absorption component in 
the column in front of the star as: 



[^wind \ 
Fc\m /„ 



Aph(v) 



)(1 



(6) 



5.1.3. The emission term 

Assuming that emission in the Ho' line comes entirely from 
spontaneous emission, thus presently neglecting any chromo- 
sphere, we can describe the emission contribution at a particular 
velocity as: 



^wind \ _ f=l^^ 
\ ^ctm /y ^ctm 



(7) 



where ;?spon is the power emitted from a shell of unit thickness 
integrated over all velocities, f^^on is the fraction of that shell at 
observed radial velocity v and Pctm is the power emitted from the 
stellar continuum (for normalisation). This latter term is approx- 
imated as a blackbody, thus Pctm is given by the familiar: 



2hv^ 



1 



exp 



X 27r X AnRldv, 



and the power p^^on is given as: 
Pspon = 47rr^ hvn^(r)Ai2, 



(8) 



(9) 



where A32 is the Einstein co-efficient for the Ha de-excitation 
(4.39 xlO^ s-i - Menzel & Pekeris 1935). 

The fraction /spon is slightly more complex, as the wind can- 
not be considered to be optically thin in all cases. We must ac- 
count for the absorption of the emitted light by the wind be- 
tween the point of emission and the observer, and it is compu- 
tationally convenient to merge this with the /spon term. We thus 
consider annuli of constant observed radial velocity (which in a 
non-rotating star also have a common depth in the wind from the 
observer's point of view), giving: 



yspon 



r — ^^pf 



-(v - Voo cos of- \ 
exp ^ , , I dv 



2w2(r) 



. Inr^ sin 9 
Anr^ 



d0, 



(10) 



where the Gaussian terms describe the fraction of atoms travel- 
ling in velocity dv, and the final term describes the fraction of 
the shell in that annulus. The lower integration limit (set to ex- 
clude flux emitted behind the star) is given as 0» = tan"'(r//?»). 
The term e''^^ describes the subsequent absorption of this emis- 
sion, as mentioned above. Defining an optical path, z, from the 
emitting region towards the observer, we can calculate r" (see 
Fig. 7): 



= J n2(r")g['o-dz, 

z=0 



(11) 



where the distance from the stellar centre, r" is given as: 

r" = ^Jiz - r' cos Of + (r' sin 6f, (12) 



where r' is the distance from the stellar centre to the emitting ma- 
terial. The function ^" selects the appropriate fraction of atoms 
in dz at that velocity, namely: 



1 



V27rw(r") 
with v" defined as: 



r") 



V - V cos y - Voo cos , 
where the angle ff' is given by: 
r" sinO" - rsinfl. 



(13) 



(14) 



(15) 



Using the above formulae, we can formulate a model where 
the only inputs are Vco, Tph, M, Tph, Vturb('" = R*), and the ra- 
dial dependences of temperature, density and (optionally) turbu- 
lence. The radial dependence of temperature can be considered 
to be constant for an isothermal wind and proportional to r^^ for 
an adiabatic, iso-kinetic wind. The radial dependence of density 
must be proportional to r^^ for a constant outflow velocity, from 
conservation of mass. /?» and Tph have been determined earlier 
from photometry and spectroscopy (Table [3]l. 

5.2. Limitations 

Our model suffers from several limitations. The model assumes 
an LTE approximation, whilst non-LTE effects may be antici- 
pated due to strong shocks and/or hot chromospheres. Attempts 
were made to model chromospheric heating in LTE using a wide 
range of temperature profiles up to T a: 7000 K, but the model 
was unable to reproduce the emission seen in many of the line 
profiles. Also, free-bound recombinations were not taken into 
account, which would alter the n - 2 and n = 3 level popu- 
lations. Although this is not significant at temperatures up to T 
~ 5500 K, it will play a role at higher temperatures where the 
plasma is ionised to a greater degree and is photo-electrically 
controlled (c.f. Thomas 1957; Dupree et al. 1984). We thus do 
not expect the model to accurately reproduce the emission seen 
in the line wings. 

Furthermore, although it contains some turbulence, our 
model assumes that the wind is isotropic. There is good evidence 
from both theory and observation (e.g. Asida & Tuchman 1995; 
Tuthill et al. 1999) that red giant winds should not be isotropic. 
Magnetic effects in the warmer stars, and pulsation and shocking 
effects in the cooler stars will mean that density enhancements 
are likely to exist in the wind as well, which may play an impor- 
tant part in shaping the observed line profile. 

Additionally, there are numerical limitations to consider: a 
limiting radius, R\im must be taken for integration and a suit- 
able step size must be chosen for dv, dr, dz, d9 and db. In our 
implementation, we use a constant step-size in all three param- 
eters, for ease of computation. In terms of R\\m, we recover the 
approximate true mass loss (assuming an inverse-square law for 
density and an isothermal, constant-velocity wind) by taking: 



^^true ~ ^ ^^computed 



r=R, r=R. 

(RUR,) - 1 



-'dr 



computed ~ 



Rlim/R* 



(16) 



where the factor | comes from the correction from a hydrogen 
mass-loss rate to a total mass-loss rate, assuming the wind is 



I. McDonald & J. Th. van Loon: Dust, pulsation, chromospheres and mass loss in red giants in Galactic GCs 



13 



75% atomic hydrogen by mass. This last value is subject to error, 
given possible helium enrichment and/or molecular hydrogen. 

We must also consider errors in the determination of the in- 
put values Tph and Varying the temperature by 100 K can al- 
ter the mass-loss rate by approximately a factor of two, due to the 
sensitive variation of the excited populations with temperature - 
the effects of a chromosphere are discussed later in this section. 
Similarly, the mass-loss rates we list are dependent on column 
densities in front of the star, thus they vary with the square of 
the stellar radius, and this then adds a ~2Q% uncertainty. 

Finally, there is also considerable degeneracy among the pa- 
rameters in the model, particularly between M and Voc, with 
a range of models able to produce similar spectra. Thus we 
conclude that our mass-loss rates are only reliable to an or- 
der of magnitude, and our wind and macroturbulence velocities 
(vturb) to within a few km s"^ Errors could be derived from 
Monte-Carlo analysis of the modelling procedure, but this is 
prohibitively computationally expensive. Using this model, we 
cannot claim to correctly determine mass-loss rates below about 
10"^ Mq yr"\ depending predominantly on the stellar temper- 
ature (high temperatures have higher populations \n n - 2 and 
3). 

5.3. The fitting metliod and results 

The observed spectrum of Ha is also affected by lines of other 
elements, notably the neutral cobalt line at 6563.403 A and the 
singly-ionised hafnium line at 6563.860 A. Thus, we took the 
best-fit Kurucz spectrum (as described in Section O, convolved 
it with a macroturbulent velocity (the value of this velocity is 
discussed later) and resampled it to the wavelengths of the model 
bins. The data were also resampled to the model wavelengths, 
re-normalised in flux using the velocity ranges between 55 and 
70 km s"' either side of the Ha rest wavelength for some stars 
(denoted a in Table|5]l, or re-normalised by setting the integrated 
flux of the model and normalised data to be equal (denoted b). A 
goodness-of-fit parameter was then calculated. To better fit the 
core of the line profiles, in which blue-shifted absorption is the 
most reliable tracer of mass loss, we decided upon a chi-squared 
statistic, with a weight proportional to the inverse square of the 
observed flux. Method b produces better fits to the Ho- profiles of 
stars with moderately-strong molecular lines, though we chose 
to accept the model with the minimum value of the chi-squared 
statistic. 

Bearing in mind the above limitations, we have fitted mod- 
els to all our targets, except those with the most extreme emis- 
sion. This was done using a downhill-simplex search in four di- 
mensions, using the parameters log(M), log(voo), log(Tph) and 
log(vturb)- The first vertex of the initial simplex had values of 
-6.5, 1, and 0.5, repectively; the other four vertices had these 
parameters changed in turn to -5.5, 1.33, -0.2 and 1.0, respec- 
tively, with the other parameters left at the original values. We 
also set dv^4km s-\ Ru^ = 4/?„ dr = 0.035/?., dz = 0.07/?,, 
d9 - nl25 and db - R,/25. Fits were also attempted using the 
following parameters, which have smaller step sizes: dv -2 km 
s-i, = 8/;., dr = O.Q2R^, dz = 0.04/?., d9 = 7r/40 and 
db - R»/25, and a larger value for but they could not re- 
produce the observations as well in the majority of cases. We 
suspect that this is due to the poorly-modelled cobalt line, which 
does not appear in most cases to be as strong as the Kurucz spec- 
trum would suggest, even when plausable forms of line broad- 
ening are taken into account. A few individual cases required 



special starting conditions in order to reproduce a good fit to the 
data, notably NGC 362 o 10, which appears to have a net inflow. 

In all cases, we chose an isothermal wind. The mass-loss 
rates required by winds with radially-decreasing temperatures 
were unreasonably high if the winds cooled by much more than 
~1000 K by 4 The macroturbulent velocity was not varied 
throughout the wind. 

Fig. [8] and Table |5] show the results of the model fitting. 
Several of the cooler stars could not be fit due to strong inter- 
ference from molecular lines, though the wind signature is often 
seen in their Ha line profiles. Additionally, NGC362o01 and 
o06, 47Tucx08, NGC 6388x03 and wCenxOl were not cor- 
rectly modelled due to their strong emission. 

In general, the absorption profiles otherwise match the ob- 
servations reasonably well. The exception being the strength of 
the other lines in the profiles. In particular, the cobalt line (at 
-1-27 km s"' with respect to Ha) is much weaker in all the ob- 
served profiles (except the M 15 stars) than the atlas9 model 
predicts. The same is true to a lesser extent of the hafnium line 
(at -t-48 km s"' with respect to Ha). Furthermore, it is interest- 
ing to note the variations in the observed line strengths when 
compared to the atlas9 model in NGC 6388. While much of 
this is likely attributable to our relatively coarse temperature grid 
(250 K spacing), it is also possible that some temperatures have 
been inaccurately determined, or that the stars have substantially 
non-solar relative abundances. Detailed modelling of the photo- 
spheric spectra with phoenix (Hauschildt et al. 1997; Baron & 
Hauschildt 1998; Hauschildt et al. 2001) will be the subject of a 
forthcoming paper. 

Two cases are worth highlighting here. All of the stars for 
which we have computed profiles show a significant amount of 
mass loss, with the exception of NGC 362 o 1 and M 1 5 xO 1 . 

NGC 362 olO cannot be modelled with an expanding wind, 
and has a Ha line core that is red-shifted beyond any error in 
radial velocity. We find it to be accreting material at a substan- 
tial rate. It is not known to pulsate strongly, nor does it exhibit 
a significantly red (J-K) colour, but it may be that it has only 
recently started to pulsate and we observe it in the contracting 
phase of its pulsation cycle where mass may have the ability to 
flow inwards. Material ejected at a velocity of ~ 10 km s"' during 
half of the pulsation cycle will reach a height of < 10 R*. The 
escape velocity of NGC 362 o 10 is one of the highest in our sam- 
ple, and will drop to < 10 km s ' only at a distance of ~76 R*. 
If the driving mechanism responsible for removing the ejected 
material from the system is too weak in this star, then the ejected 
material is bound to fall back. Indeed, pulsation-related radial 
velocity variations of this order have been seen both in field gi- 
ants (Lebzelter 1999) and globular cluster giants (Lebzelter et 
al. 2005), with wind models suggesting that inflowing material 
in the wind is also possible (Lebzelter et al. 2002). However, it 
remains a possibility that some of the blue wing of the Ha line 
is filled in by emission, moving the absorption core redwards. 
A similar phenomenon may also be present in NGC 6388 o02, if 
the emission in this object is due to a chromosphere. However, 
given the high luminosity and low temperature of this object, it 
seems more likely that the emission seen here is due to pulsation. 

Conversely, M 15 xOl appears to have both very little emis- 
sion and very little (though still clearly non-zero) mass loss. 
It perhaps represents the most quiescent of our sample, which, 
given its low metallicity, is interesting. M 15x02, on the other 
hand, displays both emission wings and a strong blue-shifted 
outflow signature, so the metallicity dependence is far from 
clear-cut. 
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Fig. 8. Results of model fitting to the Ho; lines. X-axis is velocity from -70 to +70 km s ^ y-axis is normalised flux from zero to 1.4 


times the continuum. The solid lines show the UVES spectrum; the long-dashed lines show the model fit using method (a); and the 



short-dashed lines show method (b). For NGC 362 olO, only method (a) was calculated, and this was done for an inflowing wind. 
The cross represents zero velocity and the continuum flux level. The profiles are discussed in the text. 
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Table 5. Results of fitting the Ha line profiles with the wind 
model described in Section 5. 



ID 


M 


Voo 












(Moyr-') 


(km s~ 


) 


(kms"') 












NGC 362 








o02 


7.6x10"'*: 


> 8.4 


1.93 


9.2 


240 


a 


o03 


1.5x10"^: 


> 17.4 


2.38 


7.8 


432 


a 


o04 


4.3x10-'': 


> 6.7 


1.96 


10.3 


386 


a 


o05a 


1.8x10"* 


7.2 


1.60 


13.0 


604 


b 


o05b 


1.3x10"'^ 


5.9 


1.53 


15.2 


495 


b 


o07 


5.7x10"'': 


> 12.4 


1.78 


8.0 


427 


a 


o08 


1.5x10"^ 


12.3 


1.94 


7.8 


507 


a 


o09 


1.9x10"': 


> 12.9 


1.74 


9.1 


233 


a 


olO 


-1.2x10"' 


-5.9 


1.12 


19.3 


167 


c 


xOl 


3.7x10"*: 


> 6.9 


1.60 


9.7 


660 


b 


x02a 


9.9x10"* 


12.1 


1.63 


11.4 


373 


a 


x02b 


1.3x10"' 


16.1 


0.84 


11.5 


350 


a 


x02c 


1.3x10"' 


10.7 


1.71 


14.6 


93 


a 


x03 


2.1x10"*: 


> 11.3 


1.89 


7.1 


499 


a 








47Tuc 








x02 


1.3x10"* 


14.7 


1.58 


8.2 


277 


a 


x03 


5.8x10"' 


6.5 


1.39 


7.2 


165 


a 


x04 


1.1x10* 


10.7 


1.22 


18.9 


478 


b 


x05 


2.1x10"' 


18.6 


1.90 


8.0 


938 


a 


x09 


7.2x10"* 


16.3 


1.68 


6.0 


504 


a 








NGC 6388 






x02 


3.5x10"' 


18.5 


2.04 


10.1 


412 


b 


x05 


2.4x10"' 


14.2 


1.38 


8.3 


446 


a 


o02 


1.3x10"*: 


> 1.5 


0.37 


14.3 


420 


b 


o03 


1.4x10"* 


9.5 


1.72 


8.5 


598 


a 


o05 


4.6x10"* 


10.6 


0.55 


12.6 


293 


b 


o07 


1.4x10"' 


5.3 


0.01 


34.0 


92 


a 


o08 


3.5x10"* 


7.4 


0.84 


10.7 


190 


a 


o09 


7.9x10"* 


7.0 


1.99 


12.3 


407 


a 


olO 


3.6x10"* 


10.5 


0.41 


0.7 


489 


b 


oil 


2.0x10"': 


15.4 


: 1.81 


8.6 


1066 


a 


ol2 


2.1x10"': 


17.6 


: 1.70 


11.1 


1112 


b 








M15 








xOl 


4.3x10"** 


10.0 


1.60 


14.0 


45 


a 


x02 


6.4x10"': 


> 14.9 


1.24 


13.6 


175 


a 



Notes: The fits are to an isothermal, constant velocity wind, with- 
out the inclusion of a chromosphere. The value is a relative 
goodness-of-fit parameter only. The coolest stars did not yield any 
meaningful results due to strong interference from molecular bands 
(see text). Errors are also discussed in Section [F!2l Values with a 
trailing colon are uncertain: these and lower limits to velocities are 
discussed in Section lS^l 'Method (a) is normalised using a contin- 
uum taken at ± (55 - 70) km s"' from Ha; method (b) is normalised 
using the entire line from -70 to -1-70 km s"'; method (c) uses the 
same normalisation as (a), but uses an inflowing wind (see text). 
The method used was chosen to give the best fit, i.e. lowest 



5.4. The influence of chromospheres on the spectral line 
profiles 

As noted, our simple model does not include all the physics nec- 
essary for high-energy plasmas. Thus, we cannot expect to ac- 
curately model a chromosphere and reproduce the observed line 
profile. However, we have experimented with the introduction 
of a layer in the wind with variable temperature, density and 
outflow velocity. By inserting a near-stationary layer with the 
temperature and density enhanced relative to the wind, we can 
mimic a chromosphere and use it to qualitatively explore its ef- 
fects. 



The observed emission wings seen in several profiles can- 
not be reproduced by a localised density enhancement (with the 
corresponding change in velocity). We find that, in order to re- 
produce the emission wings, there needs to be a significant frac- 
tion of gas at a temperature well above that of the photosphere 
{T > 7000/r), meaning it will be largely ionised - i.e. a chro- 
mosphere. We find it is likely restricted to radii of < 2R» for 
stars showing only a blue emission wing, as this is likely a re- 
sult of the red wing being blocked by the stellar disk; whereas 
stars showing both red and blue emission wings are likely to also 
have a large amount of material at radii > 1 .5/?*. It seems further 
likely that stars showing only a red emission wing, or P-Cygni- 
type profile, have the blue wing reabsorbed by the material in 
front of the emitting region. 

The existence of a chromosphere can mean that the true (in- 
stantaneous) mass-loss rates differ significantly from the mass- 
loss rates we have calculated above. The existence of emission 
wings creates two competing effects; firstly, the temperature en- 
hancement increases the fraction of atoms in the n - 2 level, 
leading us to overestimate the mass-loss rate; secondly, the in- 
filling of the blue-shifted absorption core leads us to underesti- 
mate both the mass-loss rate and expansion velocity. In general 
our modelling suggests that the former term dominates, and in 
the stars showing the strongest wings, the mass-loss rate may be 
overestimated by up to a factor of roughly ten and the velocity 
underestimated by up to 10 km s"' . Weaker emission wings will 
cause less of an eff'ect. The aff'ected mass-loss rates are listed as 
uncertain, and the velocities as lower limits, in Table |5] 



6. Discussion 

6.1. Spectroscopic correlations with Infrared excess 

One of our original aims was to identify differences between 
the populations with and without mid-IR excess, with the ex- 
pectation that the excesses were likely due to dust-rich winds. 
We present some statistical properties of our sample in Table |6] 
(47 Tuc has not been included for comparison due to its small 
sample size, and the strong molecular bands and effects of 
shocking in its cooler stars). 

It appears as though there is little, if any, statistical difference 
between IR-normal and IR-excessive stars whatsoever. If we 
consider the entire sample, the temperature of the IR-excessive 
stars is slightly lower, but this appears to be due to selection ef- 
fects when making up the original sample and disappears when 
we only consider the clusters where we have both excessive and 
non-excessive objects (as in Table |6]l. Interestingly, the temper- 
atures of the two sets of stars are also identical, contrary to 
the expectation the IR-excessive stars should be cooler due to 
their advanced evolution (which leads to dust production and IR- 
excess). 

One possibility for this is that the intrinsic brightness vari- 
ability of the stars, due to their pulsation, is causing some ob- 
jects to appear to have an IR excess when they do not. As the 
photometry used to assign stars with an excess is not strictly co- 
temporal, this suggests that it is normal variability, rather than 
dust, which separates our TR excessive' stars from the rest of 
the sample. As such, we would expect there to be little or no 
difference between the parameters, as is the case. However, the 
IR excesses listed in OFFR02 are of order of a magnitude, while 
the K-band amplitudes found for these stars (Table |2]) are < 0.6 
magnitudes, suggesting this cannot be the only reason. Hagen et 
al. (1983) find little correlation between mass-loss rate and dust- 
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Table 6. Properties of stars with and without IR excess in clus- 
ters NGC 362 and NGC 6388. 



Property 




Objects with: 
No excess (24) Excess (23) 


Effective 




3886 ± 280 ± 57 


3893 ± 


327 ± 


103 


lemp. {is^) 
















Ho- 








Core oil set 




-3.5 ± 4.3 ±0.9 


-0.4 ± 


9.1 ± 


2.9 


Kink offset 




1.9 ± 3.3 ±0.7 


3.5 ± 


6.1 ± 


1.9 


Vplnpitv ranc 




-5.4 ± 3.9 ±0.8 


-3.8 ± 


5.2 ± 


1.6 






Ca 11 8498 A 








Core ofiset 




1.5 ± 2.2 ±0.4 


2.7 ± 


5.3 ± 


1.7 


Kink offset 




1.9 ± 2.2 ±0.5 


0.8 ± 


1.4 ± 


0.4 


VClULiLy 




-0.4 ± 2.7 ± 0.6 


1.9 ± 


5.3 ± 


1.7 






Ca 11 8542 A 








Core offset 




-0.7 ± 2.2 ±0.5 


-0.6 ± 


5.3 ± 


1.7 


Kink offset 




2.2 ± 2.9 ±0.6 


0.0 ± 


2.3 ± 


0.7 


Vplnpitv rant 




-2.8 ± 2.6 ±0.5 


-0.6 ± 


7.3 ± 


2.3 






Ca 11 8662 A 












0.3 ± 2.3 ±0.5 


1.6 ± 


4.8 ± 


1.5 


iVillJv UllaCL 




2.0 ± 1.9 ±0.4 


1.0 ± 


2.0 ± 


0.6 


Velocity r&n^ 




-1.7 ± 2.3 ±0.5 


0.6 ± 


5.8 ± 


1.8 






Ca 11 Average 








l_Ulc UllaCL 




0.4 ± 2.0 ± 0.4 


1.3 ± 


5.0 ± 


1.6 


Kink offset 




2.0 ± 2.1 ±0.4 


0.6 ± 


1.6 ± 


0.5 


Velocity ran^ 


^e 


-1.6 ± 2.1 ±0.4 


0.6 ± 


5.9 ± 


1.9 






All objects in: 










NGC 362 (12) 


NGC 6388 (18) 


Effective 




4221 ± 186 ± 46 


3713 ± 


219 ± 


52 


lemp. yis^) 
















Ho- 








Core offset 




-3.3 ± 4.6 ±1.3 


-1.7 ± 


7.4 ± 


1.8 


Kink offset 




2.5 ± 3.6 ±1.0 


2.7 ± 


4.9 ± 


1.2 


\/plr»pit\7 v'3r\( 
VClUUlLy ld.li± 




-5.8 ± 4.5 ± 1.3 


-4.4 ± 


3.9 ± 


0.9 






Ca 11 8498 A 








Core offset 




1.9 ± 1.3 ±0.4 


2.1 ± 


4.5 ± 


1.0 


Kink offset 




3.0 ± 2.7 ±0.8 


0.9 ± 


0.9 ± 


0.2 


\/plr»piti/ rnnf 
VClUUlLy ld.li£ 




-1.1 ± 2.1 ±0.6 


1.2 ± 


4.5 ± 


1.1 






Ca 11 8542 A 








Core offset 




0.2 ± 1.9 ±0.6 


-0.5 ± 


3.8 ± 


0.9 


Kink offset 




3.4 ± 2.9 ±0.8 


0.4 ± 


2.5 ± 


0.6 


Velocity ranj 




-3.2 ± 3.2 ±0.9 


-0.9 ± 


5.2 ± 


1.2 






Ca 11 8662 A 








Core offset 




0.9 ± 1.3 ±0.4 


0.7 ± 


4.3 ± 


1.0 


Kink offset 




3.2 ± 1.8 ±0.5 


0.7 ± 


1.1 ± 


0.3 


Velocity ranj 




-2.3 ± 2.0 ±0.6 


0.0 ± 


4.6 ± 


1.1 






Ca 11 Average 








Core offset 




1.0 ± 1.2 ±0.3 


0.7 ± 


4.1 ± 


1.0 


Kink offset 




3.2 ± 2.4 ±0.7 


0.7 ± 


1.2 ± 


0.3 


Velocity ranj 




-2.2 ± 2.2 ±0.7 


0.1 ± 


4.5 ± 


1.1 



Notes: The other clusters have been excluded to reduce selection 
effects. The first error column gives the standard deviation from 
the mean of the sample, the second error is the error in the mean. 
Velocity range is defined by kink velocity (or redmost velocity 
when the kink is not present) minus core velocity. Units are km 
s"', except where noted. 



to-gas ratio, thus we should perhaps not expect any differences 
in the IR-excessive stars, which presumably contain more dust. 

It should be noted that the objects with excess show higher 
variances in velocity offsets and ranges, particularly in the core 
velocity. Considering the period-magnitude diagram (Fig.|5]l we 
find that the objects with longer periods are the objects with 
stronger pulsations (Table |2]i. All the objects for which we have 



variabihty information are (barring NGC362o01) identified as 
IR-excessive in the isocam photometry, as are the few most lu- 
minous and coolest stars (see Fig.|4|i. 

The chromospheric activity appears more pronounced in IR- 
normal stars, but this bias is most likely due to selection effects, 
since the stars in NGC 362 (which contains most stars showing 
chromospheric activity) are hotter and more luminous, therefore 
are more likely to have chromosphere-driven winds, rather than 
pulsation-driven winds. However, if IR excess is linked with dust 
formation, then it would not be surprising that dust formation is 
inhibited while the hot chromosphere is still in place. 

These factors suggest that those objects labelled IR- 
excessive stars are, by and large, more active objects with greater 
variability - hence more evolved objects. Further evidence for 
this comes from Fig. |4] which shows that the most luminous, 
coolest objects are all IR-excessive. We do not find any signif- 
icant correlation between period and temperature, nor between 
core asymmetry and luminosity, nor core asymmetry and tem- 
perature. It is important to realise that the stars are within a fac- 
tor ~4 of the RGB tip: stars lower on the RGB are likely to pul- 
sate less, and dust formation at lower luminosities may be more 
stochastic and a by-product of outflows, rather than an integral 
part of driving a wind. 

Also in Table |6] we compare the two clusters, which are 
at significantly different metallicities (Table [TJ. Given our rel- 
atively small sample size, there are no statistically significant 
variations between the clusters in the Ha line, despite the clear 
difference in stellar temperature. This suggests that temperature 
is not a strong influence on the outflow velocity. 

6.2. Line core velocities 

Fig. |9]shows the variation of hydrogen versus calcium line core 
velocities, with respect to the stellar rest velocity and the 'kink' 
velocity described previously. Clearly a correlation exists in both 
cases, although the calcium lines are only perturbed by around 
30% of the velocity of the hydrogen lines. This suggests that the 
lower layers of the atmosphere, where the calcium lines form, 
have not yet been accelerated to the terminal velocity. The two 
outliers showing redshifted Ha cores in the lower panel of Fig. |9] 
are NGC 362 oOl and NGC 6388 x03, both of which show strong 
emission in Ha. 

Though their results are mostly for the Na D fines rather than 
Ha and Ca II, Bates et al. (1993) find that low-luminosity stars 
show only small core velocity shifts, though these increase be- 
yond a luminosity of s; 1000 L© and become more scattered. 
Most of our sample is above 1000 Lq, thus we cannot verify 
this claim. We do see considerable scatter, but do not find any 
correlation between line core velocity and luminosity within the 
luminosity range of our sample. 

6.3. Line profile shapes 

Fig.|6]clearly shows unusually deep line profiles. To investigate 
this, and the other line profile features already noted, the at- 
LAs9 model corresponding to the temperature and metallicity of 
the star was subtracted from the Ha profile. Interpolations were 
made to more accurately model the correct metallicity. These re- 
sults are presented in Figs.fTOl&fTTl 

Two things are immediately apparent from these figures: 
firstly, that the absorption components are indeed much deeper 
than the atlas9 modelling would suggest; and secondly, that 
emission in both red and blue components persists. We have dis- 
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Fig. 9. Plots of Ha core velocity against calcium triplet core ve- 
locity (averaged over the three lines). Velocity is taken with re- 
spect to (a) the stellar radial velocity and (b) with respect to the 
line 'kink'; see Sect. 4 for details. The dotted lines indicate zero 
and equal velocities for both the Ha and Ca 11 lines. 



cussed a chromospheric model for producing the emission. The 
deep absorption could arguably be produced by the same phe- 
nomenon: from an optically thick chromosphere. 

Another facet of Figs. [10] & [HI is the progression of emis- 
sion and line asymmetry with spectral type and metallicity. Note 
first that in the figures, the emission spike present in most stars 
in NGC6388 (and possibly 47Tuc) at 6563.4 A co-incides with 
a Co I line. This line has been grossly over-subtracted, suggest- 
ing that there may be a cobalt deficiency in these objects, or an 
error in our temperature calculations for intermediate (~4000 K) 
temperatures. However, examination of the original data (Fig. 
|6|l suggests that at least some part of this emission may be still 
associated with the Ha line. 

Looking within clusters, blue asymmetry and line depth ap- 
pear to increase as luminosity increases towards 2400 L© (or as 
temperature decreases towards ~3900 K). As luminosity con- 
tinues to increase and temperature decrease, profiles become 
more variable. Similarly, broad emission wings only appear in 
stars below ^ 2400 Lq (above ~3900 K), and in cool stars 
is only present if they show high photometric variability and 
pulsation shocks. The emission wings are only seen clearly in 
NGC362 and M 15, the low-metallicity clusters (and possibly 
NGC6388o02, though at this luminosity and temperature this 
could be due to pulsation), with little or no emission seen in 




6560 6562 6564 

A [Al 

Fig. 10. Ha line profiles of the target stars, sorted by temper- 
ature (temperature increases upwards), with the appropriate at- 
LAs9 model spectrum subtracted. Spectra have been normalised 
to their continuum fluxes as in Figs. [T]&|6]and shifted to lab- 
oratory wavelengths. Spectra for 47Tucx06, x07 and x08, and 
wCenxOl have not been shown, as they are too cool to be ac- 
curately matched to an atlas9 model. The four panels show ob- 
jects in different clusters, with metallicity decreasing from left 
to right. A Co i line at 6563.4 A appears to have been oversub- 
tracted in NGC6388, giving the false impression of red-shifted 
emission. 



the higher-metaflicity NGC6388 and 47Tuc. While our low- 
metallicity sample are generally warmer stars where chromo- 
spheres may be expected, and where there are fewer strong spec- 
tral lines to mask their signature, the complete lack of apparent 
line wing emission in the high-metallicity stars is striking. This 
may mean that mass loss from metal-poor stars is dominated by 
their chromospheres longer in the evolution of the star Note that 
the deep Ha line, blue asymmetry and emission appear to dis- 
appear for the cool stars (which are not the most luminous) in 
NGC6388. We remind the reader that we cannot readily differ- 
entiate whether stars below »2400 L© are on the RGB or AGB, 
which could confuse these results. 



6.4. Mass-loss rates and wind velocities 

The calculated mass-loss rates can be found in Table |2l which 
also includes the following semi-empirically-determined mass- 
loss rates from Reimers (1975), SC05 and Nieuwenhuijzen & de 
Jager (1990, hereafter NdJ90), respectively (in Mq yr '): 



Mr75 = 4 X lO-'^T?— ^ 
M 



LR 



(17) 
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Table 7. Estimated mass-loss rates and equivalent widths (EW) 
of the Ha line profile emission contribution. 



6560 6562 6564 6560 6562 6564 6560 6562 6564 6560 6562 6564 
X [A] X [1] A \k] \ \k] 

Fig. 11. As Fig.[TO]but with spectra sorted by luminosity (most 
luminous at top). 



Msco5=8xlO-"*— f^^r 
" M UOOOK/ 



1 + 



8o 



4300^. 



MNdj90 = 9.631 X IO^'^l'-'^^M^ '^T?"-^', 



(18) 
(19) 



for luminosity L, radius R, mass M in solar units and stellar and 
solar gravities g, and (see Table [3] where M = 0.8 has 
been adopted); and a semi-empirically determined constant rj, 
which is usually in the range 0.4 < 77 < 3 (Stancliffe & Jeffery 
2007), with ?] ~ 0.5 being most favoured (SC05). Although val- 
ues up to // = 10 have been used effectively in specific cases 
(e.g. Straniero et al. 1997 in modelling for intermediate-mass 
field stars), whereas SC05 express t] in terms of Tgff and g, based 
on a heuristic model for the chromospheric driving of mass- 
loss. In this work, we assume rj - 1, as it provides a closer 
match to the SC05 results. We also compare our results to those 
of OFFR02, who determine mass-loss rates near-independently 
(save for identical photometry with which we have calculated 
the stellar parameters) and Gratton (1983), which is discussed 
later in this section. OFFR02 use dust masses inferred from in- 
frared colours and scale them to estimate the total mass loss, via 
a gas-to-dust ratio, for which they adopt: 



Mgas/Mdust = 200 X 10 



-([Fe/H]+0.76) 



(20) 



Figure[T2]shows a comparison between the first five methods 
in Table Q Although they weakly follow the predicted correla- 
tion with stellar parameters, all our targets clearly lie above the 
mass-loss rate predicted by all three semi-empirical relations. 
However, our mass-loss rates agree well with those determined 
by OFFR02 - the observed scatter is concordant with our ex- 
pected errors. As the OFFR02 results are derived independently 



ID 


EW 






M(10- 










(A) 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 








NGC 362 








oOl 






1.3 


1.1 


0.26 






0O2 


0.14 


0.8 


0.3 


0.15 


0.05 




3.2 


o03 


0.14 


1.5 


1.2 


1.0 


0.24 




5.8 


o04 


0.10 


4.3 


0.6 


0.34 


0.10 




3.5 


o05a 


0.11 


18 


1.0 


1.3 


0.19 




4.6 


o05b 




13 


1.4 


0.99 


0.29 






0O6 






1.4 


1.2 


0.28 






o07 


0.18 


57 


1.2 


0.98 


0.24 




6.1 


0O8 


0.13 


154 


1.1 


0.88 


0.22 




5.2 


o09 


0.11 


1.9 


0.15 


0.07 


0.02 




2.2 


olO 




-1.2 


0.13 


0.07 


0.02 






xOl 


0.10 


37 


0.8 


0.59 


0.16 


4.4 


4.1 


x02a 




99 


1.4 


0.98 


0.29 






x02b 


0.16 


133 


1.6 


1.2 


0.34 




6.4 


x02c 




1.3 


0.8 


0.48 


0.16 


2.8 




x03 


0.16 


21 


1.1 


0.85 


0.21 


3.6 


5.4 








47Tuc 








xOl 






7.2 


13.8 


1.85 


2.5 




x02 




14 


0.5 


0.29 


0.08 


0.9 




x03 




5.8 


1.2 


0.96 


0.24 


1.0 




x04 




11 


0.5 


0.31 


0.09 


2.1 




x05 




2.1 


2.1 


2.2 


0.45 


11 




x06 






8.7 


17.8 


2.23 


4.8 




x07 






9.3 


21.1 


2.55 


3.4 




xOB 






19.6 


64.5 


5.73 


46fe) 




x09 




72 


2.7 


3.6 


0.66 












NGC 6388 








oOl 






1.1 


0.71 


0.19 






0O2 


0.34 


13 


3.4 


4.4 


0.78 




20.1 


o03 




14 


0.5 


0.28 


0.08 






o04 






1.2 


0.88 


0.22 






o05 




46 


0.5 


0.30 


0.09 






0O6 






1.3 


1.1 


0.26 






o07 




1.4 


0.5 


0.31 


0.08 






0O8 




35 


1.1 


0.81 


0.21 






o09 




79 


0.4 


0.19 


0.07 






olO 




36 


0.5 


0.28 


0.09 






oil 




2.0 


0.7 


0.38 


0.12 






0I2 




2.1 


1.0 


0.62 


0.17 






xOl 






0.9 


0.52 


0.15 


1.4 




x02 




3.5 


1.0 


0.69 


0.18 


1.4 




x03 






2.1 


2.09 


0.44 


2.2 




x04 






0.8 


0.44 


0.13 


1.2 




x05 




2.4 


1.4 


1.20 


0.29 


6.8 




x06 






2.0 


1.87 


0.41 


2.3 










M15 








xOl 


0.03 


0.4 


0.09 


0.04 


0.01 


3.3 


0.9 


x02 


0.06 


6.4 


0.23 


0.12 


0.03 


4.3 


2.0 










Cen 








xOl 






2.3 


1.67 


0.38 







Equivalent widths are determined using residuals from the 
model fits. Mass-loss rates determined using or quoted from: 
(a) our model (section 5); (b) Reimers (1975), (c) Schroder 
& Cuntz (2005); (d) Nieuwenhuijzen & de Jager (1990); (e) 
OFFR02; (f) Gratton (1983), assuming a shell temperature of 
8500 K and using the peak of the emission wing as v^; (g) 
rate for 47Tucx08 from Frogel & Elias (1988), using similar 
methods to OFFR02. The errors on these values are large and 
discussed in the text. 
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Fig. 12. Comparison of estimated mass-loss rates from our modelling, Reimers (1975), OFFR02, SC05 and NdJ90. Key: NGC 362: 
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mass-loss rates are very weakly correlated with, but significantly higher than those predicted by Reimers, SC05 and NdJ90, but in 
reasonable agreement with those observationally determined by OFFR02, which are determined through measured dust masses and 
dust temperatures. 



from ours, it seems probable that our mass-loss rates - at least 
for the IR-excessive stars - are roughly correct. 

The reason behind the lack of agreement with the semi- 
empirical relations is unclear, but could be influenced by sev- 
eral factors. Though Reimers' formula focussed on giants, the 
mass-loss rates for all three semi-empirical relations are mainly 
derived from, on average, more massive, more luminous and hot- 
ter stars than we see here. The Reimers and SC05 models also 
do not make allowances for winds influenced by dust or pulsa- 
tion, which are undoubtedly important in some of our sample. It 
has also been suggested that the NdJ90 models use an unphys- 
ical prescription based on a mix of mass-loss rates from dust- 
driven and non-dust-driven sources (Schroder & Cuntz 2007). 
Furthermore, the models make no account of individual varia- 
tions between stars outside of the parameters they use, such as 
pulsation phase and composition; they are also calibrated using 
only galactic stars, and relations for stars at near-solar metallic- 
ities may not necessarily hold for low-metallicity environments 
such as these. Origlia et al.'s values, however, are calculated for 
each star individually using the dusty model (Ivezic et al. 1999). 
The systematically higher mass-loss rates that we find with re- 
spect to the predictions from semi-empirical relations cannot be 
accounted for merely by a bias in the target selection procedure 



(Section 2), and it suggests instead that the efficiency of mass 
loss for low-mass, metal-poor red giants at the tip of the RGB 
and AGB is underestimated by the semi-empirical relations. 

In exploring the correlation between our mass-loss rates and 
velocities and various other stellar parameters, we find little cor- 
relation among any of them. Perhaps the strongest correlation is 
between mass-loss rate and luminosity (Fig.[T3]l, which shows a 
general increase in mass-loss with luminosity (hence evolution), 
in a similar fashion as is seen in more massive dust-enshrouded 
AGB stars in the Magellanic Clouds (e.g. van Loon et al. 2005). 
Note that the coolest and most luminous stars which show the 
strongest emission have not been modelled, but are expected 
to have mass-loss rates near the top of this distribution, which 
would further enhance this correlation. 

We also present the relationships between modelled wind ve- 
locity and computed escape velocity (Fig.O and wind momen- 
tum versus eff'ective temperature (Fig.[T5]l after JS91. In the first 
instance, like JS91, we do not find any correlation between the 
two parameters, but we do find similiar average velocities around 
10 km s ^ Our wind momenta follow roughly the relation they 
suggest, but the momentum we derive appears to be several times 
lower for a given temperature. We suggest that this is due to the 
fact that our stars are less massive (likely ~0.5-0.8 Mq, depend- 



20 



I. McDonald & J. Th. van Loon: Dust, pulsation, chromospheres and mass loss in red giants in Galactic GCs 



10" 



10'- 



to" 



10'' 



10-' 









+ + 

+ 

+ x 

++ 


X 


- 






+ 








□ 










+ 








+ 










+ 


+ 


+ + 




B 




. . 1 







1000 

Luminosity (solar units) 



Fig. 13. Variation of mass-loss rate with luminosity. Key is as in 
Figure [T2] 



ing mostly on mass-loss history), wheras some of their sample 
could be up to 4 Mq. We do not find that wind velocity decreases 
with later spectral type (c.f. point four in JSQl's conclusions). 
We also cannot support the claim that wind velocity increases 
with luminosity, as Cahn & Wyatt (1978) suggest for field stars. 
It is also noteworthy that the velocities we find are significantly 
lower than the escape velocity of the star Given our results are 
most sensitive to the material within the first few stellar radii, this 
suggests that the material must receive an accelerating force long 
after it leaves the stellar surface. Once again, with their strong 
pulsations the coolest (unmodelled) stars are likely to have large 
wind momenta and thus contribute to the correlation observed. 
We do not find any correlation between metallicity and wind ve- 
locity, nor turbulent velocity (c.f. Gratton et al. 1984). 

Both Figs. [T3] & [T4l show broad scatter, the former particu- 
larly towards higher luminosities. This scatter is greater than that 
expected from our errors, and could possibly represent variation 
of mass-loss rate with the phase of the pulsation cycle or from 
differences between the RGB and AGB. In both figures there are 
hints of a bimodality in both NGC 362 and NGC 6388 individ- 
ually, and in the sample as a whole. A histogram in Figure [T6l 
shows that the distribution is indeed slightly bimodal, avoiding 
the mass-loss region of ~10"^ Mq yr"', though the strength of 
this bimodality depends sensitively on the size of the bins. 

The cause of this bimodality is not clear It may correspond 
to an 'on' and an 'off' phase of mass-loss as the star expands 
and contracts, or could be related to chromospheric instability. 
The lack of correlation with luminosity suggests that this is not a 
product of two different evolutionary phases, though considering 
there is only one star in the low-mass-loss region that is clearly 
above the RGB tip (several others could be below it within the er- 
rors), the bimodality may reflect a difference in mass-loss on the 
RGB and AGB. Were this the case, the Reimers and SC05 mass- 
loss laws would still hold for the weaker RGB winds, where dust 
and pulsation may not be important, though they would not ap- 
ply to the stronger, dusty and/or pulsation-driven winds of the 
AGB stars. 

The caveat to our results is that it becomes difficult to explain 
the high mass-loss rates derived for these temperatures unless all 
the stars for which we have calculated a mass-loss rate are in a 
brief phase of particularly high mass loss. Using Cioni et al.'s 
models, we find that an (initially) 0.8 M© star loses 0.13 Mq 
on the portion of the AGB above the luminosity of the RGB tip 
(ss 1900-2200 Lq, depending on metallicity) over a period of 
around 2.5 Myr, placing a limit on the average mass-loss rate of 

5.2 xlO"^ Mq yr ' during that phase. The average value we find 6.5. Shell temperatures and chromospheres 
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Fig. 14. Variation of wind velocity with escape velocity. Key is 
as in Figure[T2] Attached triangles denote lower limits. 



We find no significant deviation between our mass-loss rates and 
those of Origlia et al. among the clusters, suggesting that their 
scaling of gas-to-dust with metallicity may be correct. 



for these stars from modelling is 4.2 xlO"^ Mq yr"', and may 
be larger, depending on the mass-loss rates of the most luminous 
stars that were not modelled. We also note that the highest mass- 
loss rates we find (> 4 x 10"^ Mq yr"') are among stars without 
IR excess or observed pulsation (except for 47 Tuc x09), but also 
without strong evidence of chromospheric activity (line wing 
emission) in all but two of the less extreme cases (NGC 362 o07 
and xOl). These stars are mostly around spectra types K1-K3 
and may have derived temperatures which are too low for their 
spectral type (see Fig.|2]i; the later type stars this applies to are all 
in the metal-rich NGC 6388 and the Ha line may be significantly 
affected by molecular lines that were not correctly modelled. 

If we assume the presence of a chromosphere, we can re- 
duce our derived mass-loss rates to lie more in line with theory. 
However, this may mean that the results given by OFFR02 are 
also too high. This could be due to incorrect assumptions on 
their part, particularly in the assumed gas-to-dust ratio (although 
that would require near-solar values even in metal-poor stars). 



In order to further explore the properties of the chromosphere 
and outflow, we also compute mass-loss rates from the Ha emis- 
sion wings using the formula modified from Cohen (1976): 



M = -i]v^RJRsWAf-^e-^-^''^\ 



(21) 



(in Mo yr"') where is the radius of the shell in stellar radii 
(assumed to be around IR,, though this value is not critical when 
compared to the uncertainty in the other parameters), Wa is the 
equivalent width (in A) of the emission line (see Table |7|i, and 
T4 is the temperature of the shell in units of lO"* K. The constant 
T] has been given as 6.6 x 10"'*' by Gratton (1983), which we 
adopt here; and as 2.4 x 10"" by Cacciari & Freeman (1983) 
(both derived using G & K stars). The origin of the difference 
between these two values is unclear The factor | arises from the 
conversion from hydrogen mass-loss rate to total mass-loss rate 
(ignoring molecular hydrogen). 
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Table 8. Comparison of wind velocities and escape velocities. 
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Fig. 15. Variation of wind momentum with stellar effective tem- 
perature. Key is as in Figure [121 The dotted line represents the 
fit by Judge & Stencel (1991) for field giants. 
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Fig. 16. Histogram of mass-loss rate. The entire sample is repre- 
sented in light grey, with the corresponding Poissonian errors; 
NGC362 is shown on the left of each bin in dark grey and 
NGC6388 is shown in black on the right of each bin. Note the 
bimodality seen both in each cluster and the sample as a whole. 



It should be noted that Cohen's model (used by Gratton and 
Cacciari & Freeman) bases the outflow velocity on the peak of 
the emission wing. These velocities are typically a factor ~4 
larger than the velocities we obtain from modelling the absorp- 
tion profile. We follow their procedure and estimate the velocity 
from the emission wings, which are listed in Table [8] Gratton 
(1983) calculates T4 as: 



T4 = 10" 



1 -H(r./2.2x 104)ln[;;2/^^]' 



(22) 



produces temperatures of typically around 2400-2550 K for all 
stars due to the terms in the denominator It is difficult to see 
how a shell at this temperature could produce the substantial 
amount of emission seen. On the contrary, exploration of our 
model suggests that the shell temperatures must be > 6000 K in 
order to produce the observed line profiles, though we cannot de- 
termine this accurately using only the physics we have incorpo- 
rated into the model. This latter figure would be in better agree- 
ment with other predictions, such as those of Falceta-Gongalves 



ID 



Velocity (km s ) 





Core 


Model 


Wing 




3R. 




(a) 


(b) 


(c) 


(d) 


(e) 






NGC 362 






002 


1.9 


8.4 


43.5 


72 


41 


oVj 


— U. / 


17.4 


43 


j4 


T 1 


o04 


2.4 


6.7 


44.5 


64 


37 


o05a 


6.5 


7.2 


41 


56 


32 


o07 


6.2 


12.4 


40.5 


54 


31 


0O8 


7.9 


12.3 


41 


54 


31 


o09 


5.4 


12.9 


43 


82 


47 


xOl 


3.8 


6.9 


42.5 


58 


34 


x02b 


12.5 


16.1 


43.5 


53 


31 


x03 


3.5 


11.3 


40 


55 


32 






NGC 6388 






0O2 


-4.7 


1.5 


60 


43 


25 






M15 








xOl 


4.2 


10.0 


44.5 


92 


53 


x02 


12.9 


14.9 


47 


77 


45 



Velocities: (a) blueshifted Ha line core velocity; (b) v„ from 
our model; (c) line wing emission velocity (averaged where 
both wings show emission); (d) escape velocity at surface; (e) 
escape velocity at 3 R,. The emission in NGC 6388 o02 may be 
due to pulsation, but it is included here for completeness. 



et al. (2006). We assume a shell temperature of 8500 K (follow- 
ing Dupree et al. 1984). The mass-loss rates thus computed with 
Eq. (ISTT i are listed in Table |7] and plotted in Fig. [TT] where we 
compare them with the mass-loss rates from our model. There is 
a fair correlation and agreement between Gratton's and our mod- 
els, except that our model suggests an order of magnitude higher 
rates for some of the more extreme stars. 

The velocities we list in Table[8]are from the line wing emis- 
sion peaks, for consistency with the previous works. They are 
identical within the ~5-10% error to those velocities found from 
measuring the peak of the residual spectrum after subtraction of 
our model. Interestingly, these velocities are comparable to the 
escape velocity found at around ~3 R*. Assuming the shell is at 
roughly this radius from the star (both Gratton and Cacciari & 
Freeman assume 2 R*), this weakens the need for a mechanism 
to drive the wind beyond this radius, thus allowing a situation 
where mass-loss can occur solely through the chromosphere. 
In the case of NGC 6388 o02, pulsation may already be strong 
enough to eject matter from the surface. 

Additionally, an anti-coiTelation exists between escape ve- 
locity and wind velocity as measured by our model (Fig. [T4l l; 
though it does not exist in the sample as a whole, the anti- 
correlation may be weakly seen on a per-cluster basis. This 
would suggest that gravity inhibits the wind acceleration. 

6.6. Chromospheric disruption and tine metaiiicity 
dependence of mass ioss 

Our observations suggest that, as luminosity increases and the 
stellar surface cools, the dominance of chromospheric activity 
is replaced by radial pulsation. Unless by coincidence the chro- 
mosphere switches off at the point where significant pulsation 
begins, it seems likely that pulsation may play a more active 
role in diminishing the chromospheric control on atmospheric 
behaviour Either pulsation disrupts the mechanism responsible 
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Fig. 17. Comparison of mass-loss rates between the model of 
Gratton (1983) and our model. The dotted line indicates a one- 
to-one correlation. The expansion velocities used in the calcula- 
tion of the Gratton mass-loss rate are taken from the emission 
peak, in keeping with their analysis. 



for heating the chromosphere, or the pulsation-induced shocks 
overrun the chromosphere and displace the hot gas to further 
out in the wind where the shocks dissipate, turning the expand- 
ing column in front of the star (which provides absorption and 
hence blue asymmetry in the profile) into an emitting column of 
material. This leads to a weakening of the line core and the onset 
of the type of profiles seen in NGC 362 oOl and NGC 6388 x03, 
and then those seen in 47 Tuc x08 and u) Cen xOl . 

The transition from chromospheric activity to pulsation- 
related activity is not absolute, and the wide variety of profiles 
may represent either pulsation or a chromosphere dominating 
activity, depending on the phase of the pulsation cycle. This 
transition appears to happen around spectral type K3-K5.5 in 
NGC 362, with pulsation-related shocking becoming noticable 
at around spectral types M1-M2 in NGC 6388. This is much 
earlier than the ~M5 spectral type for field stars found by JS91, 
which is highly suggestive of a metallicity relation among the 
three. This corresponds to temperatures of ~3500 K in the field 
stars, -3700 K in NGC 6388 and -3800-4000 K in NGC 362 
(see Fig.|2]). As can be seen from Figs. [TO]&[TT] it is not clear 
from our sample whether the transition is more dependent on 
temperature or luminosity, though these temperatures would all 
roughly correspond to the RGB tip (see Fig. |4]i. It is possible 
that our lack of differentiation between the RGB and AGB could 
be influencing our conclusions, though pulsation-induced shocks 
are only visible above the RGB tip. Interestingly, these temper- 
atures also mark the divergence of the temperature vs. spectral 
type relations of Levesque et al. (Fig.|2l), though it is not obvious 
why the two phenomena should be linked. 

The fractional occurence of Wa emission was investigated 
by Gratton et al. (1984), who find that 70% of stars above 500 
Lq (with no metallicity variation) and ~80% of stars above 2000 
Lq show Wa emission in some form, though they did not dif- 
ferentiate between chromospheric-like emission and pulsation- 
related emission. Our data show 37% and 42%, respectively, 
for the same statistics. We have only one star (M 15x01) be- 
low 500 and, of the remainder, ten (22%) are presumed to 
have chromosphere-induced emission and seven (15%) exhibit 
pulsation-induced emission (on the assumption NGC 6388 o02 
exhibits pulsation-induced emission), with a strong metallicity 



bias towards the presence of emission in the two most metal- 
poor clusters. Clearly the results are disparate. This may be a re- 
sult of a higher average metallicity of our sample (Gratton et al.'s 
sample averages [Fe/H] - -1.35 whereas ours averages [Fe/H] 
= -0.91), which in turn could affect the visibility of Wa emis- 
sion through masking by stronger atomic or molecular lines due 
to higher abundances and lower temperatures. Their sample only 
contains two stars from 47 Tuc and none from NGC 362 or 6388, 
so we cannot perform a comparative study on individual clusters. 

The mass-loss rates appear to be no different for stars in 
metal-poor and metal-rich clusters with similar properties, (Fig. 
[T2] | such as luminosity (Fig. [T3] l. This correlates with the state- 
ment made by JS91, that mass-loss rate in this phase of evolution 
does not strongly depend on the driving process. 

7. Conclusions 

In this study, we have presented VLT/UVES data of a set of gi- 
ant branch stars in globular clusters, which are among the high- 
est resolution, high signal-to-noise spectra of their kind. We use 
these to characterise and quantify the outflow from their atmo- 
spheres. 

We have used Kurucz's atlas9 models to determine stel- 
lar temperature and, from this, we have calculated basic stellar 
parameters, which appear typical for red and asymptotic giant 
branch stars. 

Many of the stars we have investigated show clear emis- 
sion in Ha, most spectacularly when strong pulsation is present, 
which is seen to occur at luminosities above the RGB tip. It is 
also noted that many of the stars show strongly blue-shifted ab- 
sorption cores, suggesting bulk outflow from the stellar surface. 
This is also mirrored in some cases in the near-infrared calcium 
triplet line profiles. 

In an effort to quantify this outflow, we have constructed a 
simple model for the warmer and more metal-poor stars. We 
have calculated terminal velocities for the winds which are of 
order 10 km s ' - much lower than the escape velocity of the 
stars - and mass-loss rates of ^ 10"^ to 10"^ M© yr"', which 
lie well above theoretical expectations, but are consistent with 
the mass-loss rates derived from IR emission from circumstel- 
lar dust. These models also suggest that an emissive shell exists 
close to the stellar surface in many stars, with the wind prob- 
ably being largely isothermal beyond it. Outflow velocities of 
the emission are of order 40 km s"', and mass-loss rates derived 
from this emission are similar to the ones we obtain from mod- 
elling the absorption profile. In the most extreme cases we may 
have either over-estimated the mass-loss rate from our model, or 
under-estimated the temperature of the emitting shell. Mass-loss 
rates correlate weakly with luminosity, but stars showing strong 
IR excesses (linked with dust production) do not necessarily ex- 
hibit higher gas mass-loss rates. We find no correlation between 
mass-loss rate and metallicity. 

We suggest that the emission and mass-loss in early-type 
(^ K3-K5) giants are dominated by a warm chromospheric re- 
gion, as suggested by some studies, though late-type (^ M3) gi- 
ants have mass loss dominated by pulsation. It seems likely that 
the spectral type of the transition between the two regimes is 
metallicity-dependent, occurring at later spectral types for higher 
metallicities. The outward velocity of the warm emissive shells 
associated with the chromospheres is similar to the escape ve- 
locity at their anticipated radius (2 to 3 R*), allowing the chro- 
mosphere to be the sole driver of mass loss in these stars. 
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